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Summary of Conference Arrangements 



The First International Conference on Smarandache Notions in Number Theory was held in 
Craiova , Romania, 21- 22 August 1997. The Organizing Committee had spared no effort in 
preparing programme, lodging and conference facilities. The Conference was opened by the la 
te professor Constantin Dumitrescu*, chairman of the Organizing Committee and the initiator 
of the conference and a leading personality in Number Theory research. He welcomed all 
participants. Unfortunately professor Dumitrescu ’s state of health did not permit him to 
actively lead the conference, although he delivered his first paper later in the day and was 
present during most sessions. He requested the author of these lines to chair the first day of the 
conference, a task for which I was elected to continue for the rest of the conference. 

In view of the above it is appropriate that I express mine and the other participants gratitude to 
the organizers and in particular to the Dumitrescu family who assisted throughout with social 
and arrangements and the facilities required for the smooth running of the conference. I would 
like to pay special tribute to professor Dumitrescu’s son Antoniu Dumitrescu who presented his 
father’s second paper on his behalf. 

Unfortunately not all those who intended to participate in the conference were able to come. 
Their contributions which were submitted in advance have been gratefully received and are 
included in these proceedings. A list of participants is found on page 

A pre-conference session was held with professor V. Seleacu the day before the conference. 
This was held in french with Mrs Dumitrescu as interpreter. Prof. Seleacu showed some 
interesting work being conducted by the research group at Craiova University. Mrs Dumitrescu 
also acted actively during the conference to bridge language difficulties. 

Special thanks were expressed at the conference to Dr. F. Luca, USA, who helped during 
sessions when translation from the romanian language to english was needed. In this context 
thanks are also due to my wife Anne-Marie Rochard-Ibstedt who made my participation 
possible by helping me drive from Sweden to Paris and then across Europe to Craiova. She was 
also active during the conference in taking photos and distributing documents. 

Although united through the international language of Mathematics it was not always possible 
to penetrate presentations in such detail that extended discussions could take place after each 
session. Informal contacts between participants proved important and opportunities for this was 
given during breaks and joint dinners. 

In the concluding remarks the chairman thanked the organizers and in particular professor 
Dumitrescu for having very successfully organized this conference. It was noted that the 
presentations were not made as an end in itself but as sources for further thought and research 
in this particular area of Number Theory, n.b. the very large number of open problems and 
notions formulated by Florentin Smarandache. The hope was expressed that the conference had 
linked together researchers for continu in g exchange of views with our modem means of 
communication such as electronic mail and high speed personal computers. 

Professor Dumitrescu thanked the chairman for his work. 



Paris 26 March 1998. 



1 1949-1997, Obituary in Vol. 8 of the Smarandache Notions Journal. 




On Smarandache’s Periodic Sequences 



Preamble: 



Henry Ibstedt 



Ladies and gentlemen. 

It is for me a great honour and a great pleasure to be here at this conference to present some 
of the thoughts I have given to a few of the ideas and research suggestions given by Florentin 
Smarandache. In both of my presentations we will look at some integer sequences defined by 
Smarandache. As part of my work on this I have prepared an inventory of Smarandache sequences, 
which is probably not complete, but nevertheless it contains 133 sequences. I welcome contributions to 
complete this inventory, in which an attempt is also made to classify the sequences according to 
certain main types. 

Before giving my first presentation I would like to say a few words about what eventually 
brought me here. 

When I was young my interest in Mathematics began when I saw the beauty of Euclidean 
geometry - the rigor of a mathematical structure built on a few axioms which seemed the only ones 
that could exist That was long before I heard of the Russian mathematician Lobachevsky and 
hyperbolic geometry. But my facination for Mathematics and numbers was awoken and who can 
dispute the incredible beauty of a formula like 

e“+l=0 

and man y others. But there was also the disturbing feet that many important truths can not be 
expressed in closed formulas and that more often than not we have to resort to approximations and 
descriptions. For a long timg I was fascinated by classical mechanics. Newton’s laws provided an ideal 
framework for a great number of interesting problems. But Einstein’s theory of relativity and 
Heisenberg’s uncer tain ty relation put a stop to living and thinking in such a narrow world. Eventually 
I ended up doing computer applications in Atomic Physics. But also my geographical world became 
too narrow and I started working in developing countries in Africa, the far East and the Caribbean, far 
away from co mp uters, libraries and contact with current research. This is when I returned to numbers 
and Number Theory. In 1979, when micro computers had just started making an impact, I bought one 
and brought it with me to the depths of Africa. Since then Computer Analysis in Number Theory has 
remained my major intellectual interest and st im ula nt . 

With these words I would now like to proceed to the subject of this session. 
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On Smarandache’s Periodic Sequences 

Henry Ibstedt 



Abstract: 

This paper is based on an article in Mathematical Spectrum, Vol. 29, No 1. It concerns what happens 
when an operation applied to an n-digit integer results in an n digit integer. Since the number of n- 
digit integers is finite a repetition must occur after applying the operation a finite number of times. It 
was assumed in the above article that this would lead to a periodic sequence which is not always true 
because the process may lead to an invariant. The second problem with the initial article is that, say, 7 
is considered as 07 or 007 as the case may be in order make its reverse to be 70 or 700. However, the 
reverse of 7 is 7. In order not to loose the beauty of these sequences the author has introduced 
stringent definitions to prevent the sequences from collapse when the reversal process is carried out 

Four different operations on n-digit integers is considered. 

The Smarandache n-digit periodic sequence. Definition: Let N k be an integer of at most n digits and 
let R* be its reverse. N k 7 is defined through 

N k = R k 

The element of the sequence through 

Nic+i= I N k -Nk’ | 

where the sequence is initiated by an arbitrary n-digit integer Nj in the domain 10 n <Ni<10 n+1 . 

The Smarandache Subtraction Periodic Sequence : Definition: Let N k be a positive integer of at most n 
digits and let R* be its digital reverse. N k ’ is defined through 

N k =R k 

The element N k+I of the sequence through 

NfcH= | Nk’-c | 

where c is a positive integer. The sequence is initiated by an arbitrary positive n-digit integer Ni. It is 
obvious from the definition that O^N^IO 1 ^ 1 , which is the range of the iterating function. 

The Smarandache Multiplication Periodic Sequence : Definition: Let Ol be a fixed integer and N 0 and 
arbitrary positive integer. Nt+i is derived from N k by multiplying each digit x of N k by c retaining 
only the last digit of the product cx to become the corresponding digit of N^i . 

The Smarandache Mixed Composition Periodic Sequence: D efinition. Let N 0 be a two-digit integer 
ai*10+ao. If ai+a<)<10 thenbi= ai+ao otherwise bi= ai+ao+1. bo=|ai-aoi . We define Ni=brl0+bo. N k +i 
is derived from N k in the same way 

Starting points for loops (periodic sequences), loop length and the number of loops of each kind has 
been calculated and displayed in tabular form in all four cases. The occurrence of invariants has also 
been included. 
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Introduction 



la Mathematical Spectrum , vol 29 No 1 [1], is an article on Smarandache’s periodic sequences which 
terminates with the statement: 

“ There will always be a periodic sequence whenever we have a function f:S->S, where S is a 
finite set of positive integers and we repeat the function f/ 

We must adjust the above statement by a counterexample before we look at this interesting set of 
sequences. Consider the following trivial function f(Xk):S->S 7 where S is an ascending set of integers 
{ai, a 2 , ... a r , ... a n ): 

f Xk-i if Xk>a, 
f(XkH X k if Xfc^ar 

l X k+1 if X k <a r 

As we can see the iteration of the function f in this case converges to an invariant a,., which we may of 
course consider as a sequence (or loop) of only one member. We will however make a distinction 
between a sequence and an invariant in this paper. 

.There is one more snag to overcome. In the Smarandache sequences 05 is considered as a two-digit 
integer. The consequence of this is that 00056 is considered as as a five digit integer while 056 is 
considered as a three-digit integer. We will abolish this ambiguity, 05 is a one-digit integer and 00200 
is a three-digit integer. 

With these two remarks in mind let’s look at these sequences. There are in all four different ones 
reported in the above mentioned article in Mathematical Spectrum. The study of the first one will be 
carried out in much detail in view of the above remarks. 



la. The Two-Digit Smarandache Periodic Sequence 

It has been assumed that the definition given below leads to a repetition according to Dirichlet’s box 
principle (or the statement made above). However, as we will see, this definition leads to a collapse of 
the sequence. 

Preliminary definition. Let N k be an integer of at most two digits and let N k ’ be its digital reverse. 
We define the element N k+ i of the sequence through 



Nw-lNfc-Nk’l 



where the sequence is initiated by an arbitrary two digit integer Nj. 

Let’s write N, in the form N,=10a+b where a and b are digits. We then have 
N 2 = | lOa+b-lOb-a I =9- 1 a-b I 

The | a-b I can only assume 10 different values 0,1,2, ... ,9. This means that N 3 is generated from only 
10 different values of N 2 . Let’s first find out which two digit integers result in I a-b I - 0 , 1 , 2 , .. and 9 
respectively. 
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I a-b | Corresponding two digit integers 



0 


11 


22 


33 


44 


55 66 


77 


88 


99 














1 


10 


12 


21 


23 


32 34 


43 


45 


54 


56 


65 


67 


76 78 


87 


89 98 


2 


13 


20 


24 


31 


35 42 


46 


53 


57 


64 


68 


75 


79 86 


97 




3 


14 


25 


30 


36 


41 47 


52 


58 


63 


69 


74 


85 


96 






4 


15 


26 


37 


40 


48 51 


59 


62 


73 


84 


95 










5 


16 


27 


38 


49 


50 61 


72 


83 


94 














6 


17 


28 


39 


60 


71 82 


93 


















7 


19 


29 


70 


81 


92 




















8 


19 


80 


91 
























9 


90 




























It is now 


easy to follow the iteration of the sequence which invariably terminates in 0, table 1. 








Table 1 


. Iteration of sequence according to the preliminary definition 






la-b| 






n 2 




N 3 




N 4 






N 5 




Ns 




N* 


0 






0 
























1 






9 




0 




















2 






18 




63 




27 






45 




9 




0 


3 






27 




45 




9 






0 










4 






36 




27 




45 






9 




0 






5 






45 




9 




0 
















6 






54 




9 




0 
















7 






63 




27 




45 






9 




0 






8 






72 




45 




9 






0 










9 






81 




63 




27 






45 




9 




0 



The termination of the sequence is preceded by the one digit element 9 whose reverse is 9. The 
following definition is therefore proposed. 

Definition of Smarandache’s two-digit periodic sequence. Let N* be an integer of at most two 
digits. N k ’ is defined through 

the reverse of N k if N* is a two digit integer 

N k ’ = 1 

Nk*10 if Nk is a one digit integer 
We define the element N^i of the sequence through 
Nk+,= I Nic-Nv’ | 

where the sequence is initiated by an arbitrary two digit integer Ni with unequal digits. 

Modifying table 1 according to the above definition results in table 2. 



Table 2. Iteration of the Smarandache two digit sequence 



ia-b| 


n 2 


n 3 


N 4 


Ns 


Ns 


N* 


n 7 


1 


9 


81 


63 


27 


45 


9 




2 


18 


63 


27 


45 


9 


81 


63 


3 


27 


45 


9 


81 


63 


27 




4 


36 


27 


45 


9 


81 


63 


27 


5 


45 


9 


81 


63 


27 


45 




6 


54 


9 


81 


63 


27 


45 


9 


7 


63 


27 


45 


9 


81 


63 




8 


72 


45 


9 


81 


63 


27 


45 


9 


81 


63 


27 


45 


9 


8 ) 





7 





Conclusion: The iteration always produces a loop of length 5 which starts on the second or the third 
term of the sequence. The period is 9, 81, 63, 27, 45 or a cyclic permutation thereof. 



lb. Smarandache’s n-digit periodic sequence. 

Let’s extend the definition of the two-digit periodic sequence in the following way. 

Definition of Smarandache’s n-digit periodic sequence. 

Let N k be an integer of at most n digits and let R* be its reverse. N k ’ is defined through 

N k = R k - lcr’-f' 08 ”"*] 

We define the element N^i of the sequence through 
Nt+,= | N k -N k ’ | 

where the sequence is initiated by an arbitrary n-digit integer Ni in the do main 10 n <N,<10 n+l . It is 
obvious from the definition that 0<N k <10 aH , which is the range of the iterating function. 



Let’s consider the cases n=3, n=4, n=5 and n=6. 
n=3. 

Domain 100<N 1 ^999. . The iteration will lead to an invariant or a loop (periodic sequence) 1 . There 
are 90 symmetric integers in the domain, 101, 111, 121, ...202, 212, ..., for which N 2 =0 (invariant). 
All other initial integers iterate into various entry points of the same periodic sequence. The number 
of numbers in the domain resulting in each entry of the loop is denoted s in table 3. 

Table 3. Smarandache 3-digit periodic sequence 



s 239 11 200 240 120 

Loop 99 891 693 297 495 

It is easy to explain the relation between this loop and the loop found for n=2. Consider 
N=ao+10ai+100a 2 . From this we have IN-N , |=99|a 2 -a 0 |=l 1 -9|a 2 -ao| which is 11 times the 
corresponding expression for n=2 and as we can see this produces a 9 as middle (or first) digit in the 
sequence for n=3. 



n=4. 

Domain lOOCKNj <9999. The largest number of iterations carried out in order to reach the first 
member of the loop is 18 and it happened for Ni=1019. The iteration process ended up in the 
invariant 0 for 182 values of Nj, 90 of these are simply the symmetric integers in the domain like 
Ni=4334, 1881, 7777, etc., the other 92 are due to symmetric integers obtained after a couple of 
iterations. Iterations of the other 8818 integers in the domain result in one of the following 4 loops or 
a cyclic permutation of one of these. The number of numbers in the domain resulting in each entry of 
the loops is denoted s in table 4. 
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This is elaborated in detail in Surfing on the Ocean of Numbers by the author, Vail Univ. Press 1997. 



Table 4. Smarandache 4-digit periodic sequences 



s 


378 


259 








Loop 


2178 


6534 








s 


324 


18 


288 


2430 


310 


Loop 


90 


810 


630 


270 


450 


s 


446 


2 


449 


333 


208 


k22£ 


909 


8181 


6363 


2727 


4545 


s 


329 


11 


290 


2432 


311 


Loop 


999 


8991 


6993 


2997 


4995 



n=5. 

Domain 10000<Ni<99999. There are 900 symmetric integers in the domain. 920 integers in the 
domain iterate into the invariant 0 due to symmetries. 



Table 5. Smarandache 5-digit periodic sequences 



s 


3780 


2590 








Loop 


21978 


65934 








s 


3240 


180 


2880 


24300 


3100 


Loop 


990 


8910 


' 6930 


2970 


4950 


s 


4469 


11 


4490 


3330 


2080 


Loop 


9009 


81081 


63063 


27027 


45045 


s 


3299 


101 


2900 


24320 


3110 


Loop 


9999 


89991 


69993 


29997 


49995 



n — 6, 

Domain 100000<N, ^999999. There are 900 symmetric integers in the domain. 12767 integers in the 
domain iterate into the invariant 0 due to symmetries. The longest sequence of iterations before 
arriving at the first loop member is 53 for N=100720. The last column in table 6 shows the number of 
integers iterating into each loop. 

Table 6. Smarandache 6<Sgit periodic sequences 



Cl 


13667 " — — — — 

0 


13667 


H 


13642 12661 

13506 65340 


24093 


ffl 


12606 \26A5 ' ' — 

219970 639934 


26271 


B9 


19107 2711 7127 123320 12446 

900 3100 6300 2700 4300 


164711 


US 


26067 fa 12269 20993 Zu9 — — 

9090 31810 63630 27270 45460 


62776 


K9 


47931 14799 42603 222941 29995 ' * — 

9990 39910 69930 29970 49950 


356269 


B 


26375 11 12375 21266 4409 

90009 310061 630063 27C027 430046 


63436 






3766 


B 


*a09 It 1350 1670 510 

99099 891891 693693 297297 496495 


5260 


t 

L 10 


191 39 2648 7292 123673 12472 

99999 399991 699993 299997 499995 


165224 


s 

1 11 


162 4 1254 972 492 111 326 466 429 ' 

10909 970021 357142 616384 131868 736263 373626 ?«77a7 494605 


4726 


j 

L 12 


«23 44 146 794 377 34 525 HO 194 494 117 144 793 327 44 430 139 179 

45449 912481 724442 441834 74329 847341 703593 308236 374417 340944 318087 442724 144438 470923 341847 404294 4,7374 





2. The Smarandache Subtraction Periodic Sequence 
Definition: 

Let N k be a positive integer of at most n digits and let R* be its digital reverse. N k ’ is defined through 



9 



















N' k =R k . 10 ' , - 1 > 8 '»^1 



We define the element N k +i of the sequence through 



Nic+i= I Nt’-c | 



where c is a positive integer. The sequence is initiated by an arbitrary positive n-digit integer Ni. It is 
obvious from the definition that O^N^IO 1 * 1 , which is the range of the iterating function. 

c=l, n=2, 10^N X ^99 

When N| is of the form Ilk or 11-k-l then the iteration process results in 0, see figure la. 

Every other member of the interval 10<Nj^99 is a entry point into one of five different cyclic periodic 
sequences. Four of these are of length 18 and one of length 9 as shown in table 7 and illustrated in 
figures lb and lc, where important features of the iteration chains are shown. 



Table 7. The subtraction periodic sequence, 10<Ni<99 



Seq. No 1 


12 


20 


1 


9 


89 


97 


78 


86 


67 


75 


56 


64 


45 


53 


34 


42 


23 


31 


Seq. No 2 


13 


30 


2 


19 


90 


8 


79 


96 


68 


85 


57 


74 


46 


63 


35 


52 


24 


41 


Seq. No 3 


14 


40 


3 


29 


91 


18 


80 


7 


69 


95 


58 


84 


47 


73 


35 


62 


25 


51 


^Seq. No 4 


15 


50 


4 


39 


92 


28 


81 


17 


70 


6 


59 


94 


48 


83 


37 


72 


26 


61 


Seq. No 5 


16 


60 


5 


49 


93 


38 


82 


27 


71 





















99 

98 

88 «- 
<- (=) 67 

77 
76 
66 
65 
55 
54 
44 
43 
33 
32 
22 
21 
11 
10 
0 

Fig. la 



37 

72 

26 

61 

15 

50 

04 <- 

<-(-1) 39 -»(+9) 

92 
28 
81 
17 
70 

— > 06 4 — 

(-1)4- 59 ->(+9) 

94 
48 
83 
37 

Fig lb 



38 

82 

27 

71 

16 

60 

05 <- 

(-!)<_ 49 -*(+9) 

93 
38 



Rg lc 



l^c^9, n=2, 100^N!^999 

A computer analysis revealed a number of interesting facts concerning the application of the iterative 
function. 

There are no periodic sequences for c=l, c=2 and c =5. All iterations result in the invariant 0 after, 
sometimes, a large number of iterations. n 




For the other values of c there are always some values of Ni which do not produce periodic sequences 
but terminate on 0 instead. Those values of Nj which produce periodic sequences will either have Ni 
as the first term of the sequence or one of the values f determined by l<f<c-l as first term. There are 
only eight different possible value for the length of the loops, namely 11, 22, 33, 50, 100, 167, 189, 
200. Table 8 shows how many of the 900 initiating integers in the interval 100<Ni<999 result in each 
type of loop or invariant 0 for each value of c. 

Table 8. Loop statistics, l=iength of loop, f=ffrst term of loop 




A few examples: 

For c=2 and Ni=202 the sequence ends in the invariant 0 after only 2 iterations: 

202 200 0 

For c=9 and Ni=208 a loop is closed after only 11 iterations: 

208 793 388 874 469 955 550 46 631 127 712 208 

For c=7 and Ni=109 we have an example of the longest loop obtained. It has 200 elements and the 
loop is closed after 286 iterations: 

109 894 491 187 774 470 67 753 350 46 633 329 916 612 209 895 591 188 874 471 
167 754 450 47 733 330 26 613 309 896 691 189 974 472 267 755 550 48 833 331 
126 614 409 89 7 791 190 84 473 367 756 650 49 933 332 226 615 509 898 891 191 
184 474 467 757 750 50 43 333 326 616 609 899 991 192 284 475 567 758 850 51 
1 43 334 426 61 7 709 900 2 193 384 476 667 759 950 52 243 335 526 618 809 901 

102 194 484 477 767 760 60 53 343 336 626 619 909 902 202 195 584 478 867 761 
1 60 54 443 337 726 620 1 9 903 302 1 96 684 479 967 762 260 55 543 338 826 62 1 
119 904 402 197 784 480 77 763 360 56 643 339 926 622 219 905 502 198 884 481 
177 764 460 57 743 340 36 623 319 906 602 199 984 482 277 765 560 58 843 341 
136 624 419 907 702 200 5 493 387 776 670 69 953 352 246 635 529 918 812 211 

105 494 487 777 770 70 63 353 346 636 629 919 912 212 205 495 587 778 870 71 
1 63 354 446 637 729 920 22 213 305 496 687 779 970 72 263 355 546 638 829 921 
122 214 405 497 787 780 80 73 363 356 646 639 929 922 222 215 505 498 887 781 
180 74 463 357 746 640 39 923 322 21 6 605 499 987 782 280 75 563 358 846 641 
139 924 422 217 705 500 2 
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3. The Smarandache Multiplication Periodic Sequence 
Definition: 

Let c>l be a fixed integer and N 0 and arbitrary positive integer. N w is derived from N k by 
multiplying each digit x of N k by c retaining only the last digit of the product cx to become the 
corresponding digit of N w . 

In this case each digit position goes through a separate development without interference with the 
surrounding digits. Let’s as an example consider the third digit of a 6-digit integer for c=3. The 
iteration of the third digit follows the schema: 

xx7yyy the third digit has been arbitrarily chosen to be 7. 

xxlyyy 

xx3yyy 

xx9yyy 

xx7yyy which closes the loop for the third digit 

Let’s now consider all the digits of a six-digit integer 237456: 

237456 

691258 

873654 

419852 

237456 which closes the loop. 

The digits 5 and 0 are invariant under this iteration. All other digits have a period of 4 for c=3. 

Conclusion: Integers whose digits are all equal to 5 are invariant under the given operation. All other 
integers iterate into a loop of length 4. 

We have seen that the iteration process for each digit for a given value of c completely determines the 
iteration process for any n-digit integer. It is therefore of interest to see these single digit iteration 
sequences: 



Table 9. One-digit multiplication sequences 



c=2 


c=3 






c= 


=4 






c=5 


1 


2 


4 8 


6 


2 


i 


3 


9 7 


1 


1 


4 


6 


4 






i 


5 


5 


2 


4 


8 6 


2 




2 


6 


8 4 


2 


2 


8 


2 








2 


0 


0 


3 


6 


2 4 


8 


6 


3 


9 


7 1 


3 


3 


2 


8 


2 






3 


5 


5 


4 


8 


6 2 


4 




4 


2 


6 8 


4 


4 


6 


4 








4 


0 


0 


5 


0 


0 






5 


5 






5 


0 


0 








5 


5 




6 


2 


4 8 


6 




6 


8 


4 2 


6 


6 


4 


6 








6 


0 


0 


7 


4 


8 6 


2 


4 


7 


1 


3 9 


7 


7 


8 


2 


8 






7 


5 


5 


8 


6 


2 4 


8 




8 


4 


2 6 


8 


8 


2 


8 








8 


0 


0 


9 


8 


6 2 


4 


8 


9 


7 


1 3 


9 


9 


6 


4 


6 






9 


5 


5 




c=6 


c=7 






c : 


=8 






c=9 


1 


6 


6 






1 


7 


9 3 


1 


1 


8 


4 


2 


6 


8 


1 


9 


1 


2 


2 








2 


4 


8 6 


2 


2 


6 


8 


4 


2 




2 


8 


2 


3 


8 


8 






3 


1 


7 9 


3 


3 


4 


2 


6 


8 


4 


3 


7 


3 


4 


4 








4 


8 


6 2 


4 


4 


2 


6 


8 


4 




4 


6 


4 


5 


0 


0 






5 


5 






5 


0 


0 








5 


5 




6 


6 








6 


2 


4 8 


6 


6 


8 


4 


2 


6 




6 


4 


6 


7 


2 


2 






7 


9 


3 1 


7 


7 


6 


8 


4 


2 


6 


7 


3 


7 


8 


8 








8 


6 


2 4 


8 


8 


4 


2 


6 


8 




8 


2 


8 


9 


4 


4 






9 


3 


1 7 


9 


9 


2 


6 


8 


4 


2 


9 


1 


9 



12 















With the help of table 9 it is now easy to characterize the iteration process for each value of c. 

Integers composed of the digit 5 result in an invariant after one iteration. Apart form this we have for; 
c=2. Four term loops starting on the first or second term. 
c=3. Four term loops starting with the first term. 

c=4. Two term loops starting on the first or second term (could be called a switch or pendulum). 

c=5. Invariant after one iteration. 

c=6. Invariant after one iteration. 

c=7. Four term loop starting with the first term. 

c=8. Four term loop starting with the second term. 

c=9. Two term loops starting with the first term (pe ndul um) 



4. The Smarandache Mixed Composition Periodic Sequence 



Definition. Let N 0 be a two-digit integer a r 10+ao. If a,+ao<10 then bi= a,+ao otherwise b,= a,+a<>+l. 
b 0 =ja I -a o | . We define Ni=bi-10+bo. N w is derived from N k in the samp way. 2 

There are no invariants in this case. 36, 90, 93 and 99 produce two-element loops. The longest loops 
have 18 elements. A complete list of these periodic sequences is presented below. 

1011 20 22 40 44 80 88 70 77 50 55 10 

11 20 22 40 44 80 88 70 77 50 55 10 1 1 

12 31 42 62 84 34 71 86 52 73 14 53 82 16 75 32 51 64 12 

13 42 62 84 34 71 86 52 73 14 53 82 16 75 32 51 64 12 31 42 

14 53 82 16 75 32 51 64 12 31 42 62 84 34 71 86 52 73 14 

15 64 12 31 42 62 84 34 71 86 52 73 14 53 82 16 75 32 51 64 

16 75 32 51 64 12 31 42 62 84 34 71 86 52 73 14 53 82 16 

17 86 52 73 14 53 82 16 75 32 51 64 12 31 42 62 84 34 71 86 

18 97 72 95 54 91 18 

19 18 97 72 95 54 91 18 

202240448088 70 77 5055 10 11 20 

21 31 42 62 84 34 71 86 52 73 14 53 82 16 75 32 51 64 12 31 

22 .40 44 80 88 70 77 50 55 10 1 1 20 22 

23 51 64 12 31 42 62 84 34 71 86 52 73 14 53 82 16 75 32 51 

24 62 84 34 71 86 52 73 14 53 82 16 75 32 51 64 12 31 42 62 

25 73 14 53 82 16 75 32 51 64 12 31 42 62 84 34 71 86 52 73 

26 84 34 71 86 52 73 14 53 82 16 75 32 51 64 12 31 42 62 84 

27 95 54 91 18 97 72 95 

28 16 75 32 51 64 12 31 42 62 84 34 71 86 52 73 14 53 82 16 

29 27 95 54 91 18 97 72 95 

30 33 60 66 30 

31 42 62 84 34 71 86 52 73 14 53 82 16 75 32 51 64 12 31 

32 51 64 12 31 42 62 84 34 71 86 52 73 14 53 82 16 75 32 

33 60 66 30 33 

34 71 86 52 73 14 53 82 16 75 32 51 64 12 31 42 62 84 34 

35 82 16 75 32 51 64 12 31 42 62 84 34 71 86 52 73 14 53 82 

36 93 36 

37 14 53 82 16 75 32 51 64 12 31 42 62 84 34 71 86 52 73 14 

38 25 73 14 53 82 16 75 32 51 64 12 31 42 62 84 34 71 86 52 73 

39 36 93 36 



Formulation conveyed to the author “Let N be a two-digit number. Add the digits, and add them again if the sum is greater than 1 0. 
Also take the absolute value of their difference. These are the first and second digits of N[." 




40 44 80 88 70 77 50 55 10 11 20 22 40 

41 53 82 16 75 32 51 64 12 31 42 62 84 34 71 86 52 73 14 53 

42 62 84 34 71 86 52 73 14 53 82 16 75 32 51 64 12 31 42 

43 71 86 52 73 14 53 82 16 75 32 51 64 12 31 42 62 84 34 71 

44 80 88 70 77 50 55 10 11 20 22 40 44 

45 91 18 97 72 95 54 91 

46 12 31 42 62 84 34 71 86 52 73 14 53 82 16 75 32 51 64 12 

47 23 51 64 12 31 42 62 84 34 71 86 52 73 14 53 82 16 75 32 51 

48 34 71 86 52 73 14 53 82 16 75 32 51 64 12 31 42 62 84 34 

49 45 91 18 97 72 95 54 91 

50 55 10 1 1 20 22 40 44 80 88 70 77 50 

51 64 12 31 42 62 84 34 71 86 52 73 14 53 82 16 75 32 51 

52 73 14 53 82 16 75 32 51 64 12 31 42 62 84 34 71 86 52 

53 82 16 75 32 51 64 12 31 42 62 84 34 71 86 52 73 14 53 

54 91 18 97 72 95 54 

55 1011 20 22 40 44 80 88 70 77 50 55 

56 21 31 42 62 84 34 71 86 52 73 14 53 82 16 75 32 51 64 12 31 

57 32 51 64 12 31 42 62 84 34 71 86 52 73 14 53 82 16 75 32 

58 43 71 86 52 73 14 53 82 16 75 32 51 64 12 31 42 62 84 34 71 

59 54 91 18 97 72 95 54 

60 66 30 33 60 

61 75 32 51 64 12 31 42 62 84 34 71 86 52 73 14 53 82 16 75 

62 84 34 71 86 52 73 14 53 82 16 75 32 51 64 12 31 42 62 

63 93 36 93 

64 12 31 42 62 84 34 71 86 52 73 14 53 82 16 75 32 51 64 

65 21 31 42 62 84 34 71 86 52 73 14 53 82 16 75 32 51 64 12 31 

66 30 33 60 66 

67 41 53 82 16 75 32 51 64 12 31 42 62 84 34 71 86 52 73 14 53 

68 52 73 14 53 82 16 75 32 51 64 12 31 42 62 84 34 71 86 52 

69 63 93 36 93 

70 77 50 55 10 1 1 20 22 40 44 80 88 70 

71 86 52 73 14 53 82 16 75 32 51 64 12 31 42 62 84 34 71 

72 95 54 91 18 97 72 

73 14 53 82 16 75 32 51 64 12 31 42 62 84 34 71 86 52 73 

74 23 51 64 12 31 42 62 84 34 71 86 52 73 14 53 82 16 75 32 51 

75 32 51 64 12 31 42 62 84 34 71 86 52 73 14 53 82 16 75 

76 41 53 82 16 75 32 51 64 12 31 42 62 84 34 71 86 52 73 14 53 

77 50 55 10 11 20 22 40 44 80 88 70 77 

78 61 75 32 51 64 12 31 42 62 84 34 71 86 52 73 14 53 82 16 75 

79 72 95 54 91 18 97 72 

80 88 70 77 50 55 10 11 20 22 40 44 80 

81 97 72 95 54 91 18 97 

82 16 75 32 51 64 12 31 42 62 84 34 71 86 52 73 14 53 82 

83 25 73 14 53 82 16 75 32 51 64 12 31 42 62 84 34 71 86 52 73 

84 34 71 86 52 73 14 53 82 16 75 32 51 64 12 31 42 62 84 

85 43 71 86 52 73 14 53 82 16 75 32 51 64 12 31 42 62 84 34 71 

86 52 73 14 53 82 16 75 32 51 64 12 31 42 62 84 34 71 86 

87 61 75 32 51 64 12 31 42 62 84 34 71 86 52 73 14 53 82 16 75 

88 70 77 50 55 10 11 20 22 40 44 80 88 

89 81 97 72 95 54 91 18 97 

90 99 90 

91 18 97 72 95 54 91 

92 27 95 54 91 18 97 72 95 

93 36 93 

94 45 91 18 97 72 95 54 91 

95 54 91 18 97 72 95 

96 63 93 36 93 

97 72 95 54 91 18 97 

98 81 97 72 95 54 91 18 97 

99 90 99 




Products of Factorials in Smarandache Type Expressions 



Florian Luca 



Introduction 



In [3] and [5] the authors ask how many primes are of the form x y +y x , 
where gcd ( x , y) = 1 and x, y > 2. Moreover, Jose Castillo (see [2]) 
asks how many primes are of the Smarandache form x X2 + x 2 X3 + ... + x n Xl , 
where n > 1, x\, x 2 , ..., x n > 1 and gcd (xi, x 2 , ..., x n ) = 1 (see [9]). 

In this article we announce a lower bound for the size of the largest 
prime divisor of an expression of the type ax y + by x , where ab ^ 0, x, y > 2 
and gcd (x, y) = 1. 

For any finite extension F of Q let d F = [F : Q]. For any algebraic 
number C € F let N F (0 denote the norm of (. 

For any rational integer n let P(n) be the largest prime number P 
dividing n with the convention that P(0) = P(±l) = 1. 

Theorem 1. Let a and (3 be algebraic integers with a ■ j3 ^ 0. Let 
K = Q[a, /?]. For any two positive integers x and y let X = max (x, y). 
There exist computable positive numbers C\ and C 2 depending only on a 
and (3 such that 



p(jV K (ox y + Py x ))> Ci 




1/(<*K + 1) 



whenever x, y > 2, gcd (x, y) = 1, and X > C 2 . 

The proof of Theorem 1 uses lower bounds for linear forms in logarithms 
of algebraic numbers (see [1] and [7]) as well as an idea of Stewart (see [10]). 

Erdos and Oblath (see [4]) found all the solutions of the equation n! = 
x p ±y p with gcd (x, y) = 1 and p > 2. Moreover, the author (see [6]) showed 
that in every non-degenerate binary recurrence sequence (u n )n> o there are 
only finitely many terms which are products of factorials. 

We use Theorem 1 to show that for any two given integers a and b with 
ab ^ 0, there exist only finitely many numbers of the type ax y + by x , where 
x, y > 2 and gcd (x, y) = 1, which are products of factorials. 

Let VT be the set of all positive integers which can be written as 
products of factorials; that is 



k 

VP = {w | w = nij!, for some m 3 - > 1}. 
j= i 




Theorem 2. Let fi, ..., f 3 G Z[X, y] 6e s > 1 homogeneous polynomi- 
als of positive degrees. Assume that fi( 0, Y) ■ fi{X, 0)^0 fori = 1, s. 
Then, the equation 

/i(xf , y?) ■ ... • f 3 {xf, yf’) € VT, (1) 

with gcd (xi, yi) = 1 and x f , y £ > 2, /or i = 1, ..., s, has finitely many 
solutions xi, yi, ..., x 3 , y 3 . Moreover, there exists a computable positive 
number C depending only on the polynomials fi, ..., f 3 such that all solu- 
tions of equation (1) satisfy max (xi, yi, — , x 3 , y 3 ) < C. 

We also have the following inhomogeneous variant of theorem 2. 

Theorem 3. Let fi, ..., f 3 € Z,[X] be s > 1 polynomials of positive 
degrees. Assume that fi(0) = 1 (mod 2) fori = 1, s. Let ai, ..., a 3 and 
bi, ..., b 3 be 2s odd integers. Then, the equation 

/ 1 (a 1 xf + biyi 1 ) - ... • f 3 (a 3 xf + b s yf’) G VT, (2) 

with gcd (Xi, yi) = 1 and x u y { > 2, for i = 1, ..., s, has finitely many 
solutions Xi, yi, ..., x 3 , y 3 . Moreover, there exists a computable posi- 
tive number C depending only on the polynomials fi, ..., f 3 and the 2s 
numbers ai, b\, ..., a 3 , b 3 , such that all solutions of equation (2) satisfy 
max (xi, yi, ..., x 3 , y 3 ) <C. 

We conclude with the following computational results: 

Theorem 4. All solutions of the equation 

x y ± y x G VF with gcd (x, y) = 1 and x, y > 2, 

satisfy max (x, y) < exp 177 . 

Theorem 5. All solutions of the equation 

x y + y z + z 1 = n! with gcd (x, y, z) = 1 and x, y, z> 2, 

satisfy max (x, y, z) < exp 518. 



2. Preliminary Results 



The proofs of theorems 1-5 use estimations of linear forms in logarithms 
of algebraic numbers. 

Suppose that Ci, 0 are algebraic numbers, not 0 or 1, of heights not 
exceeding Ai, ..., Ai, respectively. We assume A m > e* for m = 1, ..., I- 
Put Q = logAi-.logAi. Let F = Q[Ci, 01- Let n u ..., m be integers, 
not all 0, and let B > max \n m \. We assume B>e 2 . The following result 
is due to Baker and Wiistholz. 
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Theorem BW ([!]). // <"*... C + 1, then 



lCr-C -II > iexp(-(16(Z + l)d F ) 2(Z+3) filogB). (3) 



In fact, Baker and Wiirtholz showed that if log £4, log£i are any- 
fixed values of the logarithms, and A = n x log £1 + — + n { log(i £ 0, then 

log | A| > -(16W f ) 2(z+2) Q log B. (4) 

Now (4) follows easily from (3) via an argument similar to the one used by 
Shorey et al. in their paper [8]. 

We also need the following p-adic analogue of theorem BW which is due 
to van der Poorten. 

Theorem vdP ([7]). Let it be a prime ideal of F lying above a prime 
integer p. Then, 

ord.(Cr-C" - 1) < (16C+ l)*) 12<,+I) ^n(logB) 2 . (5) 



The following estimations are useful in what follows. 



Lemma 1. Let n>2 be an integer, and let p <n be a prime number. 
Then 



(i) 



n n/2 < n! < n n . 



( 6 ) 



(ii) 



n 



4 (p - 1) 



< ordpn! < 



n 



P-1 



( 7 ) 



Proof. See [6]. 

Lemma 2. (1) Let s > 1 be a positive integer. Let C and X be two 
positive numbers such that C > exps and X > 1. Let y > 0 be such that 
y <C log 5 X. Then, ylogy < (C log C) log 3+1 X. 

(2) Let s > 1 be a positive integer, and let C > exp(s(s + !))• If X is 
a positive number such that X < Clog 5 X, then X < Clog 5+1 C. 

Proof. (1) Clearly, 

plogy < C log 5 X (log C + s log log X). 

It suffices to show that 

log C + s log log X < log C log X. 
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The above inequality is equivalent to 



log C(log X — 1) > slog log X. 



This last inequality is obviously satisfied since logC > s and logX > 
log log X + 1, for all X > 1. 

(2) Suppose that X > Clog s+1 C. Since s > 1 and C > exp(s(s + 1)), 

it follows that C\og 3+1 C > C > exp s. The function , is increasing 

log y 

for y > exp s. Hence, since X > C log 3+1 C, we conclude that 

Clog’+‘C , * 

log'(Clog’ +1 C) _ log'X 



The above inequality is equivalent to 



log 



3+1 , 



^log C + (s + 1) log log c) 



< 1 , 



or 



logC < 



(l + {s + l) 



log log C \ 
logC J 



S 



By taking logarithms in this last inequality we obtain 



/ lo 0- 10 0 * C \ 

loglogC<slog(l + (s + l) j <s(s + l) 



log logC 
logC 



This last inequality is equivalent to logC < s(s -1- 1), which contradicts the 
fact that C > exp(s(s + 1)). 



3- The Proofs 

The Proof of Theorem 1. ByCi, Ci, we shall denote computable 
positive numbers depending only on the numbers a and (3. Let d = die- Let 

N K (ax*+l3y x ) = p s l '-...-pi* 

where 2 < p\ < pi < ... < Pk are prime numbers. For p. = 1, d, let 
a (M) x y + be a conjugate, in K, of ax y + (3y x . Fix i = 1, ... , k. Let 

7 r be a prime ideal of K lying above Pi. We use theorem vdP to bound 
ord*. + P^y x ). We distinguish two cases: 

CASE 1. Pi | xy. Suppose, for example, that Pi | y. Since (x, y) = 1, 
it follows that pi / x. Hence, by theorem vdP, 

ord,r (a (M) x v + P (tt) y x ) = ord ff (a (M) x y ) + ord* < 
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<C X +C 2 T^~ log 4 *. (8) 

log Pi 

where C\ = d • log 2 Nk(o), and C 2 can be computed in terms of a and P 
using theorem vdP. 

CASE 2. pi / xy. In this case 

ord X (a w x y + p M y x ) =ord „{a M x y ) -hord^l - (-~j) - < 

<Ci+C 2 t^~ log 4 X. (9) 

log Pi 

Combining Case 1 and Case 2 we conclude that 

ord w (a^x y + P^y x ) < C 3 -^- log 4 X, (10) 

v 7 logpi 

where C 3 = 2 • max ( C\ , C 2 ). Hence, 

Si = ord Pi (Nk (ocx y + py *)) < C 4 ^- log 4 X. (11) 

where C 4 = dC 3 . Denote pk by P. Since pi < P for t = 1, , k, it follows, 

by formula (11), that 

fc 

log^iVK (ocx y + Py x )) < Si ■ log Pi < kC 4 P d log 4 X. (12) 

^ i= 1 



Clearly k < n(P), where ir(P) is the number of primes less than or equal to 
P. Combining inequality (12) with the prime number theorem we conclude 

that pd+l 

' d ~ (13) 



log 



^N K (ax y + py 1 )') < Cs-^—p log 4 X. 



We now use theorem BW to find a lower bound for \og(^N K (ax y + py x ) ). 
Suppose that X = y. For p = 1, ..., d, we have 



log^|a^£ y 



+ P M V X l) = log(|a (M) x y |) + log 




> 



> C& + X log 2 — C- log 3 X. 

where C 6 = min (log|a (M) | | p = 1, .... d), and C 7 can be computed using 
theorem BW. Hence, 



log (n k {ocx y + py x ) )>dC 6 +dX log 2 - dC 7 log 3 X. ( 14) 
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Let Ca = d.Ce, Cg = d log 2 , and Cio = dCj. Let also Cn be the smallest 
positive number such that 

\ Cgy > Cio log 3 y-Ca , for y > C n ■ 

Combining inequalities (13) and (14) it follows that 

pd+l 0 1 

C 5 f—=r log 4 X>Ca + C 9 X- C l0 log 3 X > -C 9 X, (15) 

log/" - 

for X > Cn. Inequality (15) clearly shows that 

p>Cn {wx) ’ 



The Proof of Theorem 2. ByCi, C 2 , we shall denote computable 
positive numbers depending only on the polynomials / 1 , f 3 - We may 

assume that fi, ..., f a are linear forms with algebraic coefficients. Let 
fi(X, Y) = caX + p{Y where (XiPi £ 0, and let K = Q[ai, Pi, ..., a,, P a ]- 
Let (xi, yi, ..., x a , y a ) be a solution of (1). Equation ( 1 ) implies that 

JJ W K (aixf + Piy f*) = nil - ... • rifc! (16) 

»=i 

We may assume that 2 < n\ < n 2 < ... < ti *,. Let X = max (x^, yi | i = 
1 , ..., s). It follows easily, by inequality (10), that 

ord 2 (n iV K (a t xf + < C t log 4 X. (17) 



Hence, 



k 

^ordsrii! <C x log 4 X. 

i=l 



By lemma 1 , it follows that 



njt < 4Ci log 4 X. (18) 

On the other hand, by theorem 1 , there exists computable constants Cgi 
and Czi, such that 




\/{d K +i) 



( 19 ) 
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whenever x £ , yi > 2, gcd (x £ , y £ ) = 1 and Xi = max (x £ , yi) > Czi ■ Let 
Ci = min (C 2 i | i = 1, s) and let Cz = max ( Czi \ i = 1, s). Suppose 

that X > Cz- From inequality (19) we conclude that 



(aixf ‘ + &yf* ) 



V=i 




i/(d K +i) 



( 20 ) 



k 

Since P | jQrii!, it follows that P < Uk- Combining inequalities (18) and 

i=i 

(20) we conclude that 



Cz 




l/(dK+l) 



<4 Cl log^X. 



( 21 ) 



Inequality (21) clearly shows that X < Ci. 

The Proof of Theorem 3. By Ci, C 2 , ..., we shall denote computable 
positive numbers depending only on the polynomials fi, ..., f s and on the 
numbers ai, fq, a s , b s . Let (x\, yi, ..., x s , y 3 ) be a solution of (2). 
Let Xi = max (x £ , yi), and let X = max (X { | i = 1, s). Finally, let 

d, 

fi(Z) = a Y[{Z-Cij). 

3 = 1 



Let K = Q[Ci,j] i<i<« , and let d = [K : Q], D = Vdi, and c = Y[d- 

i= 1 i=l 

Let 7 r be a prime ideal of K lying above 2. Let Z £ = a £ xf + hy^. We 
first bound ord* fi(Zi). First, notice that oTd K {a i b i ) = 0. Moreover, since 
/i(0) = 1 (mod 2), it follows that orcU-(Cij) = 0, for all j = 1 , ..., di. We 
distinguish 2 cases: 

CASE 1. Assume that 2 / x t y £ . Then f t {Z t ) = /,(0) = 1 (mod 2). 
Hence, ord x /i(2i) = 0. 

CASE 2. Assume that 2 | x x . In this case, ord x (y) = 0. Fix j = 
1, ..., di- Then, 

ord x (Zi - Ci j) = ord x (a,xf + (6 t yf‘ - j) ) • ( 22 ) 

Since ord x (6 i yf i ) = ord x (Ci,i) = 0, it follows, by theorem vdP, that 

ord x (6 i yf i - Qj) = ord* (btf' (C.,j) _1 - l) < Ci log 3 X { . (23) 

We distinguish 2 cases: 
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CASE 2.1. yi > Ci log 3 Xi. In this case, from formula (22) and in- 
equality (23), it follows that 

ord ff (Zi - Ci j) = ord ff (6 i t/f i - Qj) < Ci log 3 Xi . (24) 

CASE 2.2. yi < Ci log 3 Xi. In this case, 

ord ff (Zi - Ci j) = °rd x (fey? 4 + (a-ixf - Cy) ) . (25) 

Let A = dixf - Ci,i- Let H{ A) be the height of A. Clearly, 

F(A) < C 2 xf Vi . 

Hence, 

\og(H (A)) < log C 2 + dxVi log Xi < Cz + Ci log 4 Xi, 

where Cz = log C 2 , and Ci — Ci • max (di \ i = 1, ..., s). Since ordx(fe) = 
ord x(y t Xi ) = 0, it follows, by theorem vdP, that 

ordx (Zi - Cij) = ordx(l - b^y^A) < C 5 logy, log(tf(A)) log 2 x { < 

< C 5 log 3 ^ (C 3 + C 4 log 4 Xi ) . (26) 

Let Cq = 2C 4 C 5 . Also, let 

Cr = exp((C' 3 /C 4 ) 1/4 ). 

From inequalities (23) and (26), it follows that 

ordx(Z i -C ij ))<C 6 log 7 X, for X > C 7 . (27) 

Hence, 

ord 2 (n fi(Zi)) < C s log 7 X, for X > C 7 , (28) 

1=1 

where Cs = 2max (sDC 6 , c). Suppose now that 

ri/^)=nv, (29) 

i=i j = i 

where 2 < m < n 2 < ... < n*- From inequality (28) and lemma 1, it follows 
that 

jfc 

£>; < C 9 log 7 X, 
i=i 
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where C9 = 4 Cs- Hence, 

k k k k k 

log(n%!) = J2 lo ° n ^ < E n J lo S n i < (E n j) lo °(E n i) < 

j—i i=i i=i i=i i=i 

< C 9 log 7 X (log Cg + 7 log log X) , for X > C~. ( 30 ) 

Let C10 be the smallest positive number > C 7 such that 

y > log Cg + 7 log log y, for y > Cio- 

From inequality ( 30 ), it follows that 

k 

MEM < Cg log 8 X, whenever X > Cio- ( 31 ) 

i=i 

3 

We now bound log(J| fi(Zi)). Fix i = 1, s. Suppose that yi = Xi. By 

i=l 

Theorem BW, 

ailxf ) + log 

> Cn + Xi log 2 — C12 log 3 Xi, ( 32 ) 

where Cu = min (|a { | | i = 1, ..., s ), and C12 can be computed using theo- 
rem BW. Let C13 = (log 2 )/ 2 , and let C14 be the smallest positive number 
> Cio such that 

Cu + y log 2 - C12 log 3 y > C 13 y, for y > C 14 . 

From inequality ( 32 ) it follows that 

max (log \Zi\) > C13X, for X > C 14 - ( 33 ) 

On the other hand, for each i = 1, ..., s, there exists two computable 
constants Ci and C[ such that 

|/i(Zi)| > CilZil^, whenever \Z { \ > C[. 

Let C15 = min (Ci | t = 1, s ), and let C 16 = max (C/ | i = 1, ..., s). 

Finally, let Cn = max (C14, (logCi6)/Ci3). Suppose that X > Cu- Since 
|/i(Zi)| > 1 , for all i = 1 , ..., s, it follows, by inequality ( 33 ), that 

3 

logdl/^)) ^ max (i°s!/x(^)l * = i» •••» 5 ) > 

i=l 



log \Zi\= log \dixf + by? | = log (| 
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> log C15 4- max (log |Z { | | i = 1, ..., s) > log Ci 5 4- C13X , for X > C 17. 

(34) 

From equation (29) and inequalities (31) and (34), it follows that 

log Ci5 + C n x < Cg log 8 X, for X > CiT. (35) 

Inequality (35) clearly shows that X <Ci$. 

The Proof of Theorem 4. Let X = max ( x , y ). Notice that if 
x y ±y x € VT, than xy is odd. Hence, by theorem vdP, 

ord2 (x y ± y x ) = ord2(l — (Ty) x x~ y ) < 48 36 • • log 4 X. (36) 

Suppose that 

x y ±y x = ni! • ... • njfc!, (37) 

where 2 < ni < ... < nk. From inequality (36) and lemma 1 it follows that 

k k o 

£ • m < 4(53 ord 2( n i0) < 48 36 • ^ - log 4 X < 12 - 48 36 • log 4 X. (38) 

i=i i=i 

It follows, by lemma 2 (1), that 

k k k 

log(x y ± y x ) = log Tii! = 5^ logrii! < 53 ^ Io S n ' < 

i=i »=i i=i 

k k 

< (53 «i) log (53 n i) < 12 • 4836 lo g( 12 • 48 36 ) • log 5 X < 1703 • 48 36 log 5 X. 

i= 1 i = 1 

(39) 

Suppose now that X = y. Then, by theorem BW, 

log \x y ±y x | > log |x y ~y x | = log(x y ) + log |1 - y I x _y | > 

> X log 3 - log 2 - 48 10 log 3 X. (40) 

Combining inequalities (39) and (40), we conclude that 

X < X log 3 < log 2 + 48 10 log 3 X + 1703 - 48 36 log 5 X < 1704 - 48 36 log 5 X. 

(41) 

Let C = 1704 • 48 36 , and let s = 5. Since logC = log 1704 + 36 log 48 > 30, 
it follows, by lemma 2 (2), that 

X < C ■ log 6 C < 1704 - 48 36 • 147 6 . (42) 



Hence, logX < 177. 
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The Proof of Theorem 5. Suppose that (x, y, z, n ) is a solution 
of x y +y z + z x = n\, with gcd (x, y, z) = 1 and min ( x , y, z) > 1. Let 
X = max (x, y, z). We assume that logX > 519. Clearly, not all three 
numbers x, y, z can be odd. We may assume that 2 | x. In this case, both 
y and z are odd. By theorem vdP, 

ord 2 (y z + z x ) = ord 2 (l - (~y)~ z z x ) < 48 36 ^ log 4 X < 3 • 48 36 log 4 X. 

(43) 

We distinguish two cases: 

CASE 1 . y > 3 • 48 36 log 4 X. In this case, by lemma 1 , 
n/4 < ordon! = ord 2 (x y + y z + z x ) = ord 2 (y z + z x ) < 3 • 48 36 log 4 X. (44) 
Hence, 

n < 12 - 48 36 log 4 X. (45) 

By lemma 2 (1), it follows that 

n log n < 12 • 48 36 log(12 - 48 36 ) log 5 X < 1703 • 48 36 log 5 X. (46) 
We conclude that 

X log 2 < log(x y + y z + z x ) — logn! < nlogn < 1703 • 48 36 log 5 X. 

Let C = 1703 • 48 36 /log2, and let s = 5. Since logC > 30, it follows, by 
lemma 2 (2), that 

X < C log 6 C < 2457 • 48 36 • 148 6 . 

Hence, logX < 178, which is a contradiction. 

CASE 2. y < 3 • 48 36 log 4 X. Let p be a prime 'number such that p | y. 
We first show that p /x. Indeed, assume that p | x. Since gcd (x, y, z) - 1, 
it follows that p / z. QWe conclude that p / n!, therefore n < p. Hence, 

n<p<y < 3-48 36 log 4 X. 

In particular, n satisfies inequality (45). From Case 1 we know that logX < 
178, which is a contradiction. 

Suppose now that p / x. Then, by theorem vdP, 

ord p (x y + z x ) = ord p (l - (-x)~ y z x ) < 48 36 -^ log 4 X < 

< 48 36 y log 4 X < 3 - 48 72 log 8 X. (47) 

We distinguish 2 cases: 
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CASE 2.1. z > 3 • 48 72 log 8 X. In this case, by lemma 2 (1) and in- 
equality (47), 

4( p _T) < ord P n! = ord p(y" + ( xV + z x )) = 

= ord p (x y + z x )< 3 - 48 72 log 8 X. 

Hence, 

n < 12(p - 1) • 48^ log 8 X < 12y 4S 72 log 8 X < 36- 48 108 log 12 X. (48) 
From lemma 2 (1) we conclude that 

Xlog2 < log(x y + y z + z x ) = logn! < nlogn < 

< 36 • 48 108 log(36 - 48 108 ) log 13 X < 317 - 48 109 log 13 X. (49) 

Let C = 317 • 48 109 / log 2, and let s = 13. Since logC > 182, it follows, by 
lemma 2 (2), that 

X < C log 11 C < 458 • 48 109 ln 14 (458 • 48 109 ) < 458 • 48 109 • 429 14 . 

Hence, logX < 513, which is a contradiction. 

CASE 2.2. z < 3 • A3 72 log 8 X. By theorem vdP, it follows that 

ord 2 (z x + (x y + y z )) = ord 2 (l - (-x y - y z )z~ x ) < 

< 48 36 log(x y + y z ) log 3 X < 3 • 48 36 log(x y + y z ) log 3 X. (50) 

log 2 

We now bound log(x y +y z ). Let y\ = 3 • 48 36 log 4 X and z\— 3 • 48 72 log 8 X. 
Since y <yi and z < z\, it follows that 

log(x y + y z ) < log(X yi +I/ 1 1 ) < log 2 + max (yi logX, z t logyi). 

Since z x logj/i > z\ > yi logX, it follows that 

log(x y + y z ) < log 2 + z x log yi . 

From lemma 2 (1) we conclude that 

log(x y + y z ) < log 2 + z\ log yi = log 2 + ■ (y t log yi) < 

< log 2 + 48 36 log 4 X • (3 • 48 36 log (3 • 48 36 )) log 5 X < 422 • 48^ log 9 X. (51) 
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From lemma 1 and inequalities (50) and (51) it follows that 

n/4 < ord 2 n! = ord 2 (z* + (x y + y z )) < 1266 • 48 108 log 12 X. 

Hence, 

n < 5064 • 48 108 log 12 X. 

By lemma 2 (1), it follows that 

X log2 < log(x y + y z + z x ) = logn! < nlogra < 

< 5064 • 48 108 • log(5064 • 48 108 ) log 13 X < 22 - 48 111 log 13 X. 

Let C - 22 • 48 ul /log2, and let s = 13. Since logC > 182, it follows, by 
lemma 2 (2), that 

X < Clog 14 C < 22 • 48 111 • 433 14 . 

Hence, logX < 518, which is the final contradiction. 
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In this paper we seek for an answer on Smarandache type question: may one 
create the theory of Magic squares 4x4 in size without using properties of some 
concrete numerical sequences? As a main result of this theoretical investigation 
we adduce the solution of the problem on decomposing the general algebraic 
formula of Magic squares 4x4 into two complete sets of structured and four- 
component analytical formulae. 



1 Introduction 

In the general case Magic squares represent by themselves numerical or 
analytical square tables, whose elements satisfy a set of definite basic and 
additional relations. The basic relations therewith assign some constant 
property for the elements located in the rows, columns and two main diagonals 
of a square table, and additional relations, assign additional characteristics for 
some other sets of its elements. 

Judging by the given general definition of Magic squares, there is no 
difficulty in understanding that, in terms of mathematics, the problem on Magic 
squares consists of the three interrelated problems 

a) elucidate the possibility of choosing such a set of elements which would 
satisfy both the basic and all the additional characteristics of the relations; 

b) determine how many Magic squares can be constructed from the chosen 
set of elements; 

c) elaborate the practical methods for constructing these Magic squares. 

It is a traditional way to solve all mentioned problems with taking into 
account concrete properties of the numerical sequences from which the Magic 




square numbers are generated. For instance, by using this way problems was 
solved on constructing different Magic squares of natural numbers' - 5 , prime 
numbers 6 - 7 , Smarandache numbers of the 1st kind 8 and so on. Smarandache 
type question 9 arises: whether a possibility exists to construct the theory of 
Magic squares without using properties of concrete numerical sequences. The 
main goal of this paper is finding an answer on this question with respect to 
problems of constructing Magic squares 4x4 in size. In particular, in this 
investigation we 

a) describe a simple way of obtaining a general algebraic formulae of Magic 
squares 4x4, required no use of algebraic methods, and explain why in the 
general case this formula does not simplify the solution of problems on 
constructing Magic square 4x4 (Sect. 2); 

b) give a description of a set of invariant transformations of Magic squares 
4x4 (Sect. 3); 

c) adduce a general algorithm, suitable for constructing Magic squares from 
an arbitrarily given set of 16 numbers (Sect. 4); 

d) discuss the problems of constructing Magic squares from the structured 
set of 16 elements (Sect. 5); 

e) solve the problem of decomposing the general algebraic formula of 
Magic squares 4x4 into a complete set of the four-component formulae (Sect. 
6 ). 

2 Constructing the general algebraic formula of a Magic square 4x4 

A table, presented in Fig. 1(2), consists of two orthogonal diagonal Latin 
’ squares, contained symbols A, B, C, D (L,) and a, b, c, d (L 2 ). Remind 10 - " that 
two Latin squares of order n are called 

a) orthogonal if being superimposed these Latin squares form a table whose 
all n 2 elements are various; 

b) diagonal if n different elements are located not only in its rows and 
columns, but also in its two main diagonals. 

It is evident that the table 1(2) is transformed in the analytical formula of a 
Magic square 4x4 when its parameter b = 0. By using Fig. 1(2) we reveal the law 
governing the numbers of any Magic square 4x4 decomposed in two orthogonal 
diagonal Latin squares. For this aim we rearrange the sets of the symbols in the 
two-component algebraic formula 1(2) so as it is shown in Fig. 1(6). Further, a 
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table 1(6) will be called additional one. Such name of the table is justified by the 
following: 

a) the table 1(6), containing the same set of elements as the table 1(2), has 
more simple structure than the formula 1(2); 

b) there exists a simple way of passing from this table to a Magic square 
4x4: really, if one considers that Fig. 1(1) represents the enumeration of the 
cells in the table 1(6), then, for passing from this table to a Magic square it will 
be sufficient to arrange numbers in the new table 4x4 in the order 
corresponding to one in the classical square 1(5) {the Magic square of natural 
numbers from 1 to 16 }. 
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Fig. 1 . Constructing the general algebraic formula of a Magic square 4x4. 



The more simple construction of the additional table in comparison with the 
formula 1(2) and the possibility of passing from the additional table to a Magic 
square suggest solving the analogous problems on constructing the 
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corresponding additional tables instead of solving the problems on constructing 
Magic squares. Further we shall always perform this replacement of one 
problem by another. 

It is easy to establish by algebraic methods 1 - 13 that the general algebraic 
formula of Magic square of order 4 contains 8 parameters. Thus it has one 
parameter less than the two-component algebraic formula, presented in 
Fig. 1(2) with b = 0. If one takes it into account, then there appears a natural 
possibility to seek a form for the general algebraic formula of a Magic square 
4x4 basing, namely, on this two-component algebraic formula. It seems 7 that 
for introducing one more parameter in the algebraic formula 1(2) one may add 
cell-wise this formula to the Magic square, shown in Fig. 1(3) (it can be easily 
counted that the Magic constant of this square equals zero}. Thus, the general 
algebraic formula of a Magic square 4x4 {see Fig. 1(4)} is obtained as a result 
of the mentioned operation. Therefore it may be written in the simple analytical 
form 



Li + Li+W. ( 1 ) 

By analysing Fig. 1(7), in which the general formula of Magic square 4x4 is 
presented as the additional table, one may conclude, that the availability of 
eight but not of seven parameters results in a substantial violation of the simple 
regularity existing for the elements of the additional table 1(6) and by this 
reason, changing the problem on constructing a Magic square 4x4 by that on 
constructing the corresponding additional table, will not result in a facilitation 
of its solution in the general case {passing from the additional table 1(7) to the 
general algebraic formula of the Magic square 4x4 1(4) one may realise by 
means of the classical square 1(5) in the way mentioned above for the additional 
table 1(6)}. 

3 A set of invariant transformations of a Magic square 4x4 

By means of rotations by 90 degrees and mappings relative to the sides one can 
obtain from any Magic square 4x4 seven more new ones {see Fig. 2, from which 
one can judge on changes of a spatial orientation of a Magic square on the basis 
of the changes in arrangement of the symbols A, B,C and £>}. Besides for n > 4 
there exist such internal transformations ( M-transformations ) of a Magic square 
nxn (permutations of its rows and columns) by which the assigned set of 
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[n/2]{(2[n/2] - 2)!!} Magic squares nxn can be obtained 7 from one square with 
regard for rotations and mappings, where the symbol all means the product of 
all natural numbers which, firstly, are not exceeding a , and, secondly, coincide 
with it in an evenness; [a] means the integer part of a. In particular, if the cells of 
any Magic square 4x4 are enumerated so, as it is shown in Fig. 2(9), and under 
M-transformations the specific permutations of the cells of the initial square are 
meant, then, in this case the all 4 possible M-transformations of a square 4x4 
can be represented in the form of four tables, depicted in Fig. 2(9 - 12). 

It is evident, that when studying Magic squares, constructed from the same 
set of elements, it is worthwhile, to consider the only squares which can not be 
obtained from each other by rotations, mappings and M-transformations. It is 
usually said about such a family of Magic squares, that it is assigned with 
regard for rotations, mappings and M-transformations. 
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Fig. 2. A set of invariant transformations of a Magic square 4x4. 
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4 The algorithm for constructing Magic squares from an arbitrary 16 numbers 

A complete set of Magic squares 4x4 from an arbitrarily given set of 16 
numbers with regard for rotations, mappings and ^/-transformations one may 
obtain by the following algorithm 7 : 

1. Calculate the sum of all 16 numbers of the given set and, having divided 
it into 4, obtain the value of the Magic constant S of the future Magic square 
4x4; 

2. Find all possible presentations of the number S in four different terms 
each of them belonging to the given set of the numbers; 

3. If the number of various partitionings is not smaller than 14, then, using 
the obtained list of partitionings, form all possible various sets of four Magic 
rows, containing jointly 16 numbers of the given set; 

4. Among the sets of four rows, of the obtained list, find such pairs of the 
sets which satisfy the following condition: each row of the set has one number 
from various rows of the other set; 

5. It is possible to construct Magic squares 4x4 from the above mentioned 
pairs, if among the earlier found Magic rows (partitionings of the number S) 
one succeeds in finding the two rows such that 

- these rows do not contain identical numbers; 

- each row contains one by one number from various rows both of the first 
and the second set of the pair. 

When constructing Magic squares 4x4 from the obtained pairs of the sets 
consisting of four rows and the sets of the pairs of the rows corresponding to 
these pairs one should bear in mind that: 

- a four-row pair of sets (see point 4) gives a set of Magic rows and columns 
for a Magic square 4x4; 

- the found pairs of the rows (see point 5) are used for forming the Magic 
square diagonals; 

- if it is necessary to seek for Magic squares with regard for rotations, 
mappings and M-transformations, then each differing pair of rows, found for 
the given pair of sets consisting of four rows, can be utilised for construction of 
only one Magic square 4x4; 

- the algorithm can be easy realised as a computer program. 
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5 Constructing Magic squares from the structured set of 16 elements 

We shall say that a Magic square of order 4 possesses the structure (contains a 
structured set of elements) if it is possible to construct from its elements the 
eight various pairs of elements with the sum equal to 1/2 of the Magic square 
constant. For obtaining the structural pattern of a Magic square, it is sufficient 
to connect by lines each pair of the elements, forming this structure, directly in 
the Magic square. The other {implicit) way of representing the structural pattern 
of a Magic square 4x4 consists of the following: having chosen 8 various 
symbols we substitute each pair of numbers, forming the Magic square, by any 
symbol. As it has been proved by analytical methods 5 , with account for 
rotations, reflections and ^-transformations none Magic squares 4x4 exist, 
which contains in its cells 8 even and 8 odd numbers and has structure patterns 
another than ones shown in the implicit form in Fig. 3(1 - 6). In reality 7 ’ 14 this 
statement is incorrect because for such Magic squares with respect of invariant 
transformations there exist 6 more new structure plots, depicted in Fig. 3(7 - 
12 ). 
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Fig. 3. A complete set of possible structural patterns in a Magic square 4x4, 
depicted in the implicit form. 



34 





Basing on Fig. 3, for all structural patterns we shall construct a complete set of 
general structural analytical formulae. Thus, in this section we shall solve the 
problem on decomposing the general algebraic formula 1(4) in the structured 
ones. 

I. Here we present a simple method suitable for constructing general 
algebraic formulae of Magic squares possessing the structural pattern 3(1 - 4). 
Besides, we point out some singularities of these four general structured 
analytical formulae. 

As it has been established in Sect. 2 the general algebraic formula of a 
Magic square 4x4 may be represented, as the sum of the two diagonal Latin 
squares, formed by capital and small Latin letters {see Fig. 1(2)}, and the Magic 
square {Fig. 1(3)}, having a zero Magic constant. It turns out 7 that general 
structured algebraic formulae, having structural patterns 3(1 - 4), can be 
obtained if the required conditions of a structuredness at the fixed structural 
pattern are written out separately for each of the 3 tables, forming the general 
algebraic formula 1(4). In particular, diagonal Latin squares 1(2) and the Magic 
square 1(3) will have structural patterns 3(1 - 4) at the following correlations 
between their parameters {for convenience, the numbers of the written systems 
of equations are chosen so that they are identical with the numbers of structural 
patterns, shown in Fig. 3, by which these equations have been derived} : 

\.A+C=B+D, l.A+B=C+D, 3 .A+D=B+C, 4 .A+D=B+C, (2) 

c = a + d. a = c + d, c — a + d, a — c + d, 

e = 0. e = 0. 5 = 0. 

Starting from the extracted system of equations (2) one can easily prove that: 

1) The cells of an algebraic formula having the structural pattern 3(1) 
contain two sequences involving elements of the following form: 

a) a t + e, a x +a, a x + a + d, a x + d 9 a x + b, a x + a + b + e , (3) 

a x + a + b + d t a x + b + d\ 

b) a 2 , a 2 + tf, a 2 + a + d, a 2 + d-e, a 2 + b, a 2 + a + b, 
a 2 + a + b + d-e, a 2 + b + d. 

One can see from a set of sequences (3) that the regularity existing between the 
symbols of an general algebraic formula, having structural pattern 3(1), is 
complicated due to the presence of the four elements containing the symbol e 
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Fig. 4. General algebraic formulae of a Magic square 4x4 
with structural patterns 3(5 - 12). 
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{as well as in the general algebraic formula of a Magic square 4x4 shown in 
Fig. 1(4)}. Consequently, the knowledge of the regularity existing between the 
elements of the general algebraic formula with structural pattern 3(1) can not be 
of help in creating a convenient and practical algorithm for constructing 
corresponding Magic squares {as well as for the general algebraic formula 

1(4)}. 

2) The general algebraic formulae having structural patterns 3(2 - 4) are 
decomposable in sums of two diagonal Latin squares {parameters b and e are 
equal to zero). Hence, there is the simple regularity for the elements of 
additional tables of general algebraic formulae with structural patterns 3(2 - 
4) and, consequently, the 

problem on constructing such Magic squares 4x4 from a given structured set of 
1 6 elements is easy to solve by means of these three formulae. 

II. Taking into account that for structural patterns 3(1 - 4) there exists a 
simple method for constructing the general algebraic formulae (see point I) we 
present in Fig. 4 a set of 8 general algebraic formulae which possess only 
structural patterns of 3(5 - 12) {the form of representing these formulae is 
chosen so that it reveals the regularity existing between their elements}. 
Analysing the analytical formulae presented in Fig. 4 we may come to the 
following conclusions: 

1) among the all above formulae, the formulae 10 and 11 have the most 
simple structure: the set, consisting of their 16 elements, is completely defined 
by the first element of the sequence a, and the value of the parameter b; 

2) the sets of the symbols, contained in the formulae 5, 6, 7, 8 and 9, may be 
represented in the form of the two identically constructed sequences consisting 
of 8 elements {the reader can himself get assured that the same holds true also 
for general algebraic formulae possessing structural patterns 3(2 - 4)}; 

3) there are two arithmetical sequences, each containing 6 terms and having 
the same progression difference in the formula 12. Thus, the complication of the 
regularity, governing the symbols forming the algebraic formula 12, is caused 
only by four of its elements {compare with the above information concerning 
the general algebraic formula possessing structural pattern 3(1)}. 

The main conclusion which may be drawn from the above written implies 
that for constructing Magic squares having the structural patterns 3(2 - 12) it is 
preferable to use the general algebraic formulae of Magic squares 4x4, 
corresponding to these structural patterns. 
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6 Four-component algebraic formulae of Magic squares 4x4 

1. Four-component algebraic formulae of the classical Magic squares 4x4. Since a 
classical (Magic) square contains in its cells 16 different natural numbers N (1 < 
N<16) then one may write 4 ’ 12 the formula for decomposing the number N in 5 
terms: 

JV = 8a + 4b + 2c + d + 1 , (4) 

where the parameters a , b , c and d can assume only two values: either 0 or 1. By 
means of (4) any classical square 4x4 may be identically decomposed in 4 tables 
(a-, b c-, ^-components) each of them containing 8 zeros and 8 units. From 
theoretical point of view 4 there exist the only three groups of Magic squares: 

1) correct squares — all the decomposition tables are by themselves Magic 
squares: they have in all the rows, columns and in the two main diagonals by 2 
zeros and 2 units. Further such decomposition tables we shall denote as A-form. 

2) regular squares — at least one of the decomposition tables differs from 
correct one by existing at least one of the components of the formula, which is 
necessarily a regular one: each of its rows and columns contains by two zeros 
and two units, but this condition being not preserved for the main diagonal. 
Further such decomposition tables we shall denote as B-form if its both main 
diagonals contain 4 or 0 zeros (units) and C-form if its the main diagonals 
contain 1 or 3 zeros (units). 

3) irregular squares — at least one of the decomposition tables differs from 
correct and regular one by existing at least one of the components of the 
formula has one row or one column where the number of the same symbols of 
one kind is distinct from two. 

As it can be proved by analytical methods 

a) by using ,4 -forms one may construct 4 - 12 the only 1 1 different algebraic 
formulae of correct Magic squares and with account for rotations and 
reflections 7 the only 7 following 

A j A 2 A$ A x A 2 A 3 A 5 , A x A 3 A s A 7 , A x A 2 A 6 A 7 , (5) 

A 2 Ay A$ Aj, A 2 A 5 A& A 7 , A 4 A$ A$ A g, 

will be different among them, where A x -A % forms are presented in Fig. 5(1 - 8); 
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b) by using B- and C-forms one may construct 4 with account for rotations, 
reflections and M-transformations the only 15 different algebraic formulae of 
regular Magic squares 



BCAA — 


B\C\A 2 A 2j 


BxC 2 A \A^y 






BCBA — 


BxCxB 2 A 2 , 


BxCxByAi, 


BxCxByA 3 , 


B\ CxBqA^ 




B\ C 2 B 2 A x , 


B\C 2 B^A a , 


B\ C 2 B$A 4 ; 




BCBB — 


B\ CxB 2 B„ 


B { CxB 2 B a , 


BxCxB}B 4> 


BxC 2 B 2 B 5, 




B\ CfBjB^ 


B\ C 2 B 5 B 6i 







( 6 ) 



where C, - C 4 and B, - B 6 forms are presented in Fig. 5(9 - 1 8). 

c) for classical squares 4x4 the complete set of four-component algebraic 
formulae consists of algebraic formulae of the only correct and regular Magic 
squares 7 {see sets of formulae (5) and (6)} . 

2. Four-component algebraic formulae of generalised Magic squares. Denote, 
first, A -components of a correct Magic square 4x4 by the symbols F u F 2 , F 3 and 
F 4 , second, the trivial Magic square, whose 16 cells are filled with units, by the 
symbol E. As it follows from point 1, any correct classical square 4x4 can be 
represented as the sum of 5 tables (the first 3 tables should be multiplied by 8, 4 
and 2): 

8F, + 4F 2 + 2 F 3 + F 4 + E. CO 



An algebraic generalisation of this notation is the expression 

aF, + $F 2 + aFj + 8 F 4 + zE, (8) 

which represents the general recording form of a Magic square 4x4 
decomposable in the sum of the 4-th ^-components. Since the numbers of a 
classical square 4x4 may be calculated from the formula (4), the formula (8) 
obviously permits to find the symbols contained in the cells of the generalised 
correct Magic square 4x4. In particular, there exist the following relations 

1— e, 5 — £+p, 9 — £+a, 13 — E+a+p, (9) 

2 — e +5, 6 — e+P+ 8, 10 — e+a+5, 14 — e +a +P +8, 
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3 — g +cj, 7 — £ +(3 +g, 11 — £ +a +a, 15 — e + a +p +cr, 

4 — e+a+5, 8 — E+p+a+5, 12 — e+a+a+5, 16 — E+a+p +a+5. 

between natural numbers from 1 to 16 and the symbols of the generalised 
correct Magic square 4x4. 

Note that the cells of the table, shown in Fig. 6(1), contain a complete set of 
the symbols of the generalised correct Magic square 4x4. These symbols are 
arranged so that the first cell of the table contains the symbol e, the second one 
contains the symbols £ + 8 and so on. Thus, the mentioned table is additional by 
the definition and permits to construct various algebraic formulae of the 
generalised correct Magic squares of the fourth order for the assigned correct 
classical squares 4x4. 

Change the form of recording the table 6(1) by introducing the new symbols 
g, h and/with the correlations g = £ +P, h = e +a,/ = £ +a +p. The new form of 
the table is shown in Fig. 6(3). The table 6(3) makes it clear that the rows of the 
initial additional table of the generalised correct Magic square 4x4 contain the 
sequences of four numbers formed by the same regularity. Let it be also noted, 
that the new table (as it may be easily verified) completely corresponds to the 
initial one only if between its parameters s, g 9 h and / the correlation £ +/ = g 
+ h is fulfilled. Thus, for constructing concrete examples of the generalised 
correct Magic squares it is necessary to continue the search for the indicated 
sequences involving four numbers until one finds among their first terms the 
two pairs of numbers having the same sum. 

For example, the generalised correct Magic square 4x4 may be formed from 

the following eight pairs of prime numbers “the twins”: 29, 31; 59, 61; 71, 73; 
101, 103; 197, 199; 227, 229; 239, 241; 269, 271. 

The additional table, shown in Fig. 6(1) one also may use for constructing 
the algebraic formulae of the generalised regular Magic squares on the basis of 
the given classical squares 4x4. However, due to the fact that the condition of 
Magicity is not fulfilled on the diagonals of regular tables (see point 1) for 
obtaining algebraic formulae of Magic squares in this case one has to assign 
additional correlations between the parameters of the additional table. For the 
set (6) of regular formulae of the Magic square 4x4 these necessary correlations 
between the parameters of the additional table 6(1) have the form, depicted in 
Fig. 6(6). 
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17 


29 


41 


53 


47 


59 


71 


83 


227 

















e 


e + d 


e + c 


e + c + d 


e + b 


e + b + d 


e + b + c 


e + b + c+ d 


e + a 


e + a + d 


e + a + c 


e + a + c+ d 


e+a+b 


e+a+b+d 


e+a+6+c 


e+a+b+c+d 



( 1 ) ( 2 ) 



83 


113 


293 


503 


41 


71 


251 


461 


281 


311 


491 


701 


239 


269 


449 


659 



1 


7 


67 


73 


37 


43 


103 


109 


157 


163 


223 


229 


193 


197 


257 


263 



e 


e+d 


e+c 


e+c+d 


g 


g+d 


g+c 


g+c+d 


h 


h+d 


h+c 


h+c+d 




f+d 


f+c 


f+c+d 



(3) 



(4) 



(5) 



Comp 

a 


>onents 

b 


;of foi 
c 


rmula 

d 


Correlations 
between the parameters 
c = 2d , 


A 


A 


B 


C 


B 


C 


A 


A 


a = 2b, 


A 


B 


c 


A 


b = 2c, 


A 


B 


B 


C 


b = c + 2d , 


B 


B 


C 


A 


a-b + lc, 


B 


B 


B 


C 


a = c + b + Id. 



( 6 ) 

Fig. 6. Examples of constructing additional tables for the generalised correct (1-5) and 
regular (2 - AABC, 4 -ABBC, 5 - BBBQ Magic squares 4x4. 

It is noteworthy, that for the given type of a four-component regular 
formula the set of the symbols, positioned in the cells of additional tables of the 
generalised regular Magic squares 4x4, does not depend upon the form of 
additional correlations between the parameters of the additional table 6(1), in 
other words, it does not depend on the position of the C-form in the a-, b-, c-, d- 
decompositions of regular formulae. One can be immediately convinced in this 
by constructing on the basis of Fig. 6(1) all six additional tables for various 
algebraic formulae of the generalised regular Magic squares 4x4. Thus, if it is 
possible to construct one additional table for algebraic formulae of the type 
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AABC or ABBC from the given set involving arbitrary 16 numbers, then it is 
also possible to construct the other additional tables of Magic squares of the 
given type, distinct from the above constructed one by the form of additional 
conditions for the parameters of the table 6(1). With regard for the above 
stated, only 3 additional tables, filled with prime numbers, for which the reader 
is referred to Fig. 6(2, 4, 5), suffice for constructing a complete family of 
different regular Magic squares 4x4. 



Fig. 7. Examples of irregular four-component algebraic formulae of Magic squares 4x4. 

In conclusion of this section we would like to draw attention that with 
regard for rotations, mappings and M-transformations there exist 7 8 1 irregular 
four-component algebraic formulae of Magic squares 4x4. For instance, 6 
formulae of such type are presented in Fig. 7 {for splitting the formulae, shown 
in Fig. 7, in four components, it suffices to retain in the formulae, at first, only 
the digits 1 and 2 (1st component), and then, only the digits 3 and 4 (2nd 
component), etc.}. Hence, the solution of the problem on decomposing the 
general algebraic formula 1(4) into a complete set of the four-component ones 
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2468 
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has following form: there are 7 formulae for correct Magic squares 4x4 {with 
account for rotations and reflections}, 15 and 81 formulae correspondingly for 
regular and irregular Magic squares 4x4 {with account for rotations, mappings 
and Af-transformations} . Thus, it is the main conclusion of this section that the 
complete set of four-component analytical formulae of Magic squares 4x4 can 
not simplify the solution of the problem on constructing Magic squares 4x4 
from an arbitrary given set of 16 numbers but it can make so for constructing 
the generalised correct and regular Magic squares 4x4. 

7 Summary 

As it have been demonstrated in this paper discussed Smarandache type 
question - whether a possibility exists to construct the theory of Magic squares 
without using properties of concrete numerical sequences - has the positive 
answer. However, to construct this theory for Magic squares 4x4 in size, the 
new type of mathematical problems was necessary to introduce. Indeed, in 
terms of algebra, any problems on constructing Magic squares without using 
properties of concrete numerical sequences may be formulated as ones on 
composing and solving the corresponding systems of algebraic equations. Thus, 
algebraic methods can be applied for 

a) constructing the algebraic formulae of Magic squares; 

b) finding the transformations translating an algebraic formula of a Magic 
square from one form into another one; 

c) elucidating the general regularities existing between the elements of 
Magic squares; 

d) finding for an algebraic formula of a Magic square, containing m freely 
chosen parameters, the equivalent set consisting of L algebraic formulae each 
containing the number of freely chosen parameters less than m. 

The new for algebra the type of mathematical problems is presented in 
points (b) - (d). It is evident that without introducing these problems the 
algebraic methods are not effective themselves. For instance, in the common 
case (see Sect. 2) the general formula of Magic square 4x4 can not simplify the 
solution of problems on constructing Magic squares 4x4 from an arbitrary 
given set of 16 numbers. In particular, even when solution of discussed 
problems is sought by means of a computer, in calculating respect it is more 
preferable for obtaining the solution to use algorithm, described in Sect. 4, than 
one, elaborated on the base of the general formula of Magic square 4x4. But by 
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means of decomposing the general algebraic formula of Magic squares 4x4 into 
complete sets of a defined type of analytical formulae one may decrease the 
common amount of freely chosen parameters in every such formula and, 
consequently, substantially simplify the regularity existing for the elements of 
every formula. In other words, for constructing Magic squares 4x4 from an 
arbitrary given set of 16 numbers there appears a peculiar possibility of using 
the set algebraic formulae with more simple structure instead of use one 
complex algebraic formula Magic square 4x4. 
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Palindromic Numbers And Iterations of the 
Pseudo-Smarandache Function 

Charles Ashbacher 
Charles Ashbacher Technologies 
Box 294, 1 19 Northwood Drive 
Hiawatha, LA 52233 
e-mail 71603.522@compuserve.com 

In his delightful little book[l] Kenichiro Kashara introduced the Pseudo-Smarandache 
function. 

Definition: For any n > 1, the value of the Pseudo-Smarandache function Z(n) is the 
smallest integer m such that n evenly divides 

m 

E*- 

Jfc=l 



And it is well-known that the sum is equivalent to . 

Having been defined only recently, many of the properties of this function remain to be 
discovered. In this short paper, we will tentatively explore the connections between Z(n) 
and a subset of the integers known as the palindromic numbers. 

Definition: A number is said to be a palindrome if it reads the same forwards and 
backwards. Examples of palindromes are 

121,34566543, 1111111111 

There are some palindromic numbers n such that Z(n) is also palindromic. For example, 
Z(909) = 404 Z(2222) =1111 

In this paper, we will not consider the trivial cases of the single digit numbers. 

A simple computer program was used to search for values of n satisfying the above 
criteria. The range of the search was, 10 < n < 10000. Of the 189 palindromic values 
of n within that range, 37, or slightly over 19%, satisfied the criteria. 

Furthermore it is sometimes possible to repeat the function again and get another 
palindrome. 

Z(909) = 404, Z(404) = 303 

and once again, a computer program was run looking for values of n within the range 




1 < n < 10,000. Of the 37 values found in the previous test, 9 or slightly over 24%, 
exhibited the above properties. 

Repeating the program again, looking for values of n such that n, Z(n), Z(Z(n)) and 
Z(Z(Z(n))) are all palindromic, we find that of the 9 found in the previous test, 2 or 
roughly 22%, satisfy the new criteria. 

Definition: Let Z*(n) = Z(Z(Z(. . .(n)))) where the Z function is executed k times. For 
notational purposes, let Z°(n) = n. 

Modifying the computer program to search for solutions for a value of n so that n and all 
iterations Z*(n) are palindromic for i = 1, 2, 3 and 4, we find that there are no solutions in 
the range 1 < n < 10,000. Given the percentages already encountered, this should not 
be a surprise. In fact, by expanding the search up through 100,000 one solution was 
found. 

Z(86868) = 17271, Z(17271) = 2222, Z(2222) = 1111, Z(llll) = 505 

Since Z(505) = 100, this is the largest such sequence in this region. 

Computer searches for larger such sequences can be more efficiently carried out by using 
only palindromic numbers for n. 

Unsolved Question: What is the largest value of m so that for some n, Z fc (n) is a 
palindrome for all k = 0, 1, 2, . . . , m? 

Unsolved Question: Do the percentages discussed previously accurately represent the 
general case? 

Of course, an affirmative answer to the second question would mean that there is no 
largest value of m in the first. 

Conjecture: There is no largest value of m such that for some n, Z fc (n) is a palindrome 
for all k = 0, 1, 2, 3, . . . , m. 

There are solid arguments in support of the truth of this conjecture. Palindromes tend to 
be divisible by palindromic numbers, so if we take n palindromic, many of the numbers 
that it divides would also be palindromic. And that palindrome is often the product of two 
numbers, one of which is a different palindrome. Numbers like the repunits, 11... Ill 
and those with only a small number of different digits, like 1001 and 505 appeared quite 
regularly in the computer search. 
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Computational Aspect of Smarandache’s Function 



Sabin Tabirca Tatiana Tabirca 

Abstract: The note presents an algorithm for the Smarandache's function computation. The 
complexity of algorithm is studied using the main properties of function. An interesting inequality is 
found giving the complexity of the function on the set {1,2,..., n}. 

1. Introduction 

In this section, the main properties of function are reviewed. The Smarandache’s 
function notion reported for the first time in [1]. The main results concerning the 
function can be found in [1], [2], 

The function S:N -> N defined by S(n) = min{£|£!:rc} is called Smarandache’s 
function. This concept is connected with the prime number concept, because using the 
prime numbers an expresion for the function is given. The important properties that 
are used in this paper, are showed bellow. 

1. For all n eN, the inequality S(n)<n is true. When n is a prime n um ber, the 
equality is obtained . 

2- If n - p\' ■ p k f •...•/?*' is the prime number decomposition then 
S(n) = max{S(p^),S(p i f),...,S(p k m -)} . (1) 

3. The inequality S(p k ) < p-k is true, if p is a prime number. (2) 

A lot of conjectures or open problems related ot the Smarandache’s function appear in 
the number theory. Several such problems have been studied using computersand 
reported in relevant literature, e.g.[3], [4], Using the computer 




2. An Algorithm for the Smarandache function 

In the following, an algorithms for computing the function S is presented The main 
idea of the algorithm is to avoid the calculations of factorial, because the values of 
nl, for ri>10 are a very big number and cannot be calculated using a computer. 

Let (x k ) kkl a sequence of integer numbers defined by x k = k\ mod n, (V)n > 0 . Using 
the definition of sequence , the following equations can be written: 

* *.=1 (3) 

• x k + 1 = (k + 1)! mod n = (k + 1 )k\ mod n = (k + l)x k mod n . (4) 

Based on the linear equation (4), S can be calculated as 
S(n) = min{£|£!:n} = min{^:|x t = 0} . 

The algorithm for S(n) performs the computation of x l ,x 2 ,...,x k until the 0 value is 
found. The PASCAL description of this algorithm is given bellow. 

function S(n:integer):integer; 
var 

k,x:integer; 
begin 
x:=l;k:=l; 
while xoO do 
begin 

x:=x*(k+l) mod n; 
k:=k+l; 
end; 
s:=k; 
end; 

An analysis for the complexity of algorithm is presented in the following. The work- 
case complexity and the average complexity are studied. 
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Theorem 2.1 

If n = p\ x ■ p k {--...-p k „ m then the complexity of the algorithm for computing S(n) is 
0(max{p l k l ,p 2 k 2 ,...,p m kj) . 

Proof 

The algorithm computes the value S(n) generating the sequence x,,x 2 ,...,x S()l) . The 

number of operation has the complexity 0(S(n)). 

Based on (1) and (2), the following inequality holds 

S(n) = max {S(p k '),S (p 2 S (/>*” ) } < vaax{p x k x ,p 2 k 2 ,...,p m k m } . (5) 

Therefore, it can be concluded that the complexity of computing S(n) is 
0(max{p x k x ,p 2 k 2 ,...,p m k m }) . 



Other aspect of complexity is given by the average operations number. Assume the 
value S(k) is generated, where k is a number between 1 and n. This value can be 
computed with S(k) operations. Therefore the process takes S(l) operations for the 
value S(l), S(2) operations for the value S(2), ...,a.s.o. The average of the numbers 



operations is 



S = ±±S(i) 



and gives an other estimation for the complexity of 



algorithm. 

In the following, a possible upper bounds for S are sought. Obviously, S verifies the 
simple inequality 



— l " i * 

s = -Ys{i)<±Ti = 

n i-i n i= , 



n + 1 1 J_ 

2 _ 2 ” + 2 



( 6 ) 



Inequality (6) can be improved using the strong inequalities for S(i). 



Lemma 2.1. 

If i>2 is an integer number the following inequalities hold 

1. S(2i) +S(2i + 1) <3 i+1. (7) 

2. S(2i-1)+S(2i)<3i-1 . (8) 

Proof 
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Assuming i>2 folow S(2i) <7. Applying this result both inequality are true. 



Based on lemma 2.1, we can found an upper bound for the 5 better than in (6). 



Theorem 2.2 

__ | n o 12 

If n>5 is a integer number then the inequality 5 = — V 5(/) <—« + — + — is true. 

nj^ 8 4 n 

Proof 

Two cases are considered to prove the inequality. 

Case 1: n—2m 

— " ^ «n 

nS = Y, 5(0 ~ Z [5(2/ _1 ) + 5(2/ )] = 5(1) + 5(2) + 5(3) + 5(4) + £ [5(2/ - 1) + 5(2/)] 

'"1 »- 1 i=3 

Based on (8) it can be derived that 

_ n/2 nl 2 n/2 o 

n5 = 9 + 2[5(2/-l) + 5(2/)]<9 + j;(3/-l) = 2 + J;(3/-l) = 2 + --(- + l)-- = 

2 2 2 2 



/=»3 



/*3 



3 n n n 3 2 1 n 

— 2 H ( hi) = — 71 -i /2 + 2 . 

22 2 2 8 4 

3 12 

Dividing by n, we obtain the inequality S<— n + — + — . 

8 4 n 

Case 2: n=2m+l 



(9) 



(/i-I)/2 



(—0/2 



/=*3 



"5 = X 5(0 = 5(1) + ^ [5(20 + 5(2/ + 1)] = 5(2) + 5(3) + 5(4) + 5(5) + j [5(2/) + 5(2/ + 1)] 

/=l 

Using (7), it is found 
(»-!)/2 



_ (n-tyl (n-l)/2 

«5 = 14 + y [5(2/) + 5(2/ + 1)] < 14 + y (3/ + 1) =* 3 + y (3/ + 1) = 3 +• 

3 (=3 ,=i 2 2 2 2 



3 n — 1 n + 1 n-1 

+ - 



_ - . , n-1 3,1 17 _ — 3 i i, 

= 3 +—(n 2 - 1) H = —rr +— «h . Thus, 5 < — n + — h . 



3 
8 * 



8 



8 



i 11 

8 2 + 8 « 



( 10 ) 



From (9) and (10), it is found 5 < — n + min{— + — } < + — + — . 

8 4 n 2 8« 8 4« 



3. Final Remarks 
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1. Based on theorem 2.2 we can say that the average operations number for 

3 1 2 

computing the Smarandache’s function is less than —n+—+—. 



3 12 11 

2. The upper bound —n + — + — improves the previous bound —n + — . 



3. The improving process can be extended using other sort of inequalities give the 
prime numbers 2 and 3. A lemma as similar as lemma 2.1 finds the upper bounds 
the sum of sixth consecutive terms of Smarandache’s function. 



4. Using the algorithm for computing the function S, the Smarandache’s function can 
be tabulated. The values S(l),...,S(n) for all n <5000 can be found. The algorithm 
should be reviewed to be able to compute the Smarandache function for the big 



numbers. 
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THE ANALYTICAL FORMULAE YIELDING SOME 
SMARANDACHE NUMBERS AND APPLICATIONS IN 
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In this paper we study the properties of some six numerical Smarandache 
sequences. As result we present a set of analytical formulae for the computation 
of numbers in these Smarandache series and for constructing Magic squares 
3x3 in size from ^-truncated Smarandache numbers. The examples of Magic 
squares 3x3 in size of six Smarandache sequences are also adduced. 



1 Introduction 

In this paper some properties of six different Smarandache sequences of the 1st 
kind 1 are investigated. In particular, as we stated, the terms of these six 
sequences may be computed by means of one general recurrent expression 

S(n) =CT ( a » 10¥(fl ” )+a » + 1)’ . ^ 

where a n — n-th number of Smarandache sequence; q>(n) and y(a*) — some 
functions; a — an operator. For each of six Smarandache sequences, 
determined by (1), we adduce (see Sect. 2 and 3) 

a) several first numbers of the sequence; 

b) the concrete form of the analytical formula (1); 

c) the analytical formula for the calculation of n-th number in the sequence; 

d) a set of analytical formulae for constructing Magic squares 3x3 in size 

from /:-truncated Smarandache numbers; 

e) a few of concrete examples of Magic squares 3x3 in size from ^-truncated 

Smarandache numbers. 




2 Analytical formulae yielding Smarandache sequences 

1. Smarandache numbers of S r series. If <p(fl) = n +1, a = 1 and y(a„) = [lg (/2 + 
1)]+1 then of (1) the following series of the numbers, denoted as S r series, is 
generated 

1, 12, 123, 1234, 12345, 123456,... (2) 

The each number of 



X* = -1 + 



Pg(*+0,5)] 

I 

7=0 



(Jfc + l-lO 7 '), 



(3) 



corresponds to each number a k of series (2), where the notation “OgOOl” means 
integer part from decimal logarithm of y. By (3) it is easy to construct the 
analytical formula for the calculation of n-th number in the S,-series: 



a„ = 10 x " I(i/10 Xi )- (4) 

i=l 

By expressions 

A°a s = 1234...(n-l)/i; A -1 a,=234...(n-l)/i; A- 2 a„ =34... («-!)«; ... (5) 

we introduce the operator A - * {the operator of k-tnmcating the number a„ = 
1234...(n-l)/i}. Since 

A°a, = 1, A -1 a 2 - 2, A- 2 a 3 =3, ... , A~" +l a„= n, ... C 6 ) 

it is evident that by the operator A~ k from the numbers of 5, -series one may 
produce the series of the natural numbers. And, vice versa, if the operator 
A +k {the operator of k-ex tending the number n) is introduced: 

A°/i = n; A +1 n = (/i— l)n; A +2 n =(n-2)(n~l)n\ ... (7) 
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then from the series of the natural numbers one may obtain the numbers of S x - 
series: 



A °l = a 1 . A* 1 2= a 2 , A +2 3 = < 2 3 , ... , A" 1 n - a n . 



( 8 ) 



It is evident that 

a) the operators A +A: and A”* are connected with each to other. Therefore 
one may simplify their arbitrary combinations by the mathematical rule of the 
action with the power expressions {for instance, A +2 A~ 7 A +3 = A +2 - 7+3 =A" 2 }. 

b) apart from operators of ^-truncating and ^-extending of numbers from 
the left {see (5) and (7)} one may introduce operators of ^-truncating and k- 
extending of numbers from the right {for instance, (A -2 12345)= 345, but 
(12345 A- 2 ) = 123}; 

c) by means of operators of ^-truncating and ^-extending of numbers from 
the right one may represent the different relations existing between the numbers 
of 5, -series {for instance, a n ~ (ci n _ x A +I ) = (a^ +l A -1 ) and so on}. 

2. Smarandache numbers of S 2 -series. If q>(/i)=/i+ 1 ; a = y — the operator of 
mirror-symmetric extending the number ^ [(/r-f-i)/2] °f ^-series from the right 

with 1 -truncating the reflected number from the left, if n is the odd number, and 
without truncating the reflected number, if n is the even number; \|/(<2„) = 
[ lg([(n+l)/2] + 1) ] + 1, then of (1) the following series of the numbers, denoted 
as S 2 -series , is generated 

1, 11, 121, 1221, 12321, 123321, 1234321,... (9) 

The analytical formula for the calculation of n - th number in the ^-series has 
the form 



a n = [n £ ] i 10 Xl_[Is ' 1 



;=i 



1=1 



( 10 ) 



where rf = 1 + X[(/r+i)/2] + X[nl 2] “&• 

3. Smarandache numbers of S 2 -series. If cp(/z) = n+ 1 ; a = y — the operator of 
mirror-symmetric extending the number a n of Sj-series from the left with 1- 
truncating the reflected number from the right; = Pgfa + 1)] + l, then of 
(1) the following series of the numbers, denoted as S^series, is generated 
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1, 212, 32123, 4321234, 543212345, 65432123456,... 00 

The analytical formula for the calculation of /z-th number in the S 3 -series has 
the form 

a n = 10 x "{Z(ilO Xf )/10 Pgil + Iz'/10 x '}. (12) 

i=2 i=l 

4. Smarandache numbers of S 4 -series. The series of the numbers 

1, 23, 456, 7891, 23456, 789123, 4567891,... (13) 

we denote as S 4 -series. It is evident that the series of the numbers (13) is 
obtained from the infinite circular chain of the numbers 

(123456789)(123456789) ... (123456789) ... (14) 

by means of the proper truncation from the left and the right. The analytical 
formula for the calculation of /z-th number in the S 4 -series has the form 

a„ = 10" "i\l + d- 9[rf/9]}/10 i+1 , d = i + n(n-l)/2. (15) 

/= 0 

5. Smarandache numbers of S r series. The series of the numbers 

1, 12, 21, 123, 231, 312, 1234, 2341, 3412, 4123, 12345,... (16) 

we denote as S 5 -series. By (3) it is easy to construct the analytical formula for 
the calculation of /z-th number in the 5 s -series: 

a„ - I(zlO^), z = [(V8/i— 7 — l)/2], ( 17 ) 

1=1 

d= y J -Xi-(Xz + 1 1)1, t = -l + n-z{z-Y)/2. 

6. Smarandache numbers of Si-series. The series of the numbers 

12, 1342, 135642, 13578642, 13579108642,... (18) 
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we denote as S 6 -series. The analytical formula for the calculation of «-th 
number in the S 6 -series has the form 



1+2[ X , /2] £ 



a„ = {10 “2* 



1 ( 2 /— 1)/ 10 
i=l 



&2/-l / 2] 



+ 12/ 10 CZ2 ' /21 }/10 [,S 2n] . 
/= 1 



(19) 



3 Magic squares 3x3 in size from ^-truncated Smarandache numbers 

1. Magic squares 3x3 in size from k-truncated numbers of S r series. By analysing 
numbers a n of S r series one can conclude that it is impossible to construct an 
arithmetical progression from any three numbers of S r se ries. Consequently 2 , 
none Magic square 3x3 in size can be constructed from these numbers. 
However, one may truncate number a n of S r series from the left or/and the 

right by means of the operator A~ k (5). Therefore there is a possibility to 
construct the Magic squares 3x3 in size from truncated numbers of S r series. In 
particular, the analytical formula for constructing such Magic squares is 
adduced in the Fig. 1(1). If in the formula 1(1) the parameters n, r, p and q 
take, for instance, the following values: 

a) n - 7, r = 14, p = 1 and q - 3, then it generates the Magic square 3x3 
shown in the Fig. 1(2); 

b) n = 4, r = 0, p = 1 and q= 3, then the numerical square 3x3, shown in 
the Fig. 1(3), is yielded — the square 1(3) is not Magic, but it can be easy 
transformed to one by means of revising three numbers marked out by the dark 
background {the revised square see in Fig. 1(3')}; 

c) n = 4, r = 7, p - 1 and q = 3, then the numerical square 3x3, shown in 
the Fig. 1(5), is yielded — the square 1(5) also is not Magic, but it can be easy 
transformed to one by means of revising just one number marked out by the 
dark background (the revised square see in Fig. 1(5')}. 

By analysing the squares, shown in the Fig. 1(3) and 1(5), it can be easy 
understood that the analytical formula 1(1) does not hold true only in such 
cases when natural numbers, being components of numbers A ~ k a n , have 
different amount of digits. To obtain the Magic square in this case, one is to 
correct the defects of the square generated by formula 1(1) {as it made, for 
instance, in Fig. 1(30 and 1(5') for squares 1(3) and 1(5)}, or to change the 
values of parameters n, r , p and/or q correspondingly. 
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A d n+r+ p+2q 


A ~ r a . 


A r ~ P qa n+r+2p+q 


*- r - 2p °^i P 


A ~r-p-q a 

rt+r+p+q 




^“n+r+q 


A -r-2p-2q 

1 u n+r+2p+2q 


^~ r ~ P ° n+r+p 



( 1 ) 



22232425262728 

17181920212223 

18192021222324 



15161718192021 

19202122232425 

23242526272829 



20212223242526 

21222324252627 

16171819202122 



( 2 ) 







30 


18 


26 


34 


22 


42 


14 



110181 


1234 


6789 


3456 


5678 


PSISfl 


4567 


mmA 


2345 





1234 


6789 


3456 


5678 




4567 


Hum 


2345 



(3) (3') (4) 



15161718 


891011 


13141516 


10111213 


12131415 


14151617 


11121314 


:^37ril9i 


9101112 



15161718 


891011 


13141516 


10111213 


12131415 


14151617 


11121314 


SHU 


9101112 



(5) 



( 5 ') 



(r + p + 2q) n + /z(/H"l)/2 
(r + 2p) + /i(/z+I)/2 

(r + q) n + n(n+\)/2 



rn + n(n+l)/2 
(r + p + <?) /z + n(n+ 1)/2 
(r + Ip + 2q) n + /i(/H-I)/2 



(r + 2p + n + /i(«+l)/2 
(r + 2^) /i + n(/ 2 +l )/2 
(r + /?) n + /i(/rH)/2 



( 6 ) 



Fig. 1. Constructing Magic squares 3x3 from /c-truncated numbers of S[ -series. 
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It should be noted that the proper replacement of numbers A ~ k a„ in 
squares 1(2), 1(3) and 1(5) by the sum of digits of natural numbers, being 
components of A ~ k a„, gives three different Magic squares 3x3. For instance, 
the Magic square, obtained by such way from square 1(3), is depicted in 
Fig. 1(4). The explanation of this curious fact can be found in Fig. 1(6), 
presenting the analytical formula of Magic square 3x3, which is obtained 

directly from the formula 1(1) by means of the mentioned way. 

2. Magic squares 3x3 in size from k-tnmcated numbers of S 2 -series. To apply 
the methods, elaborated in point 1, for constructing Magic squares 3x3 from 
numbers of S 2 -series {see (9)}, we divide a set of S 2 -series numbers into two 

different subsequences: 

1) a,=l, 02=121, o 3 =12321, o 4 =1234321, ... 

2) £>,=11, b 2 = 1221, #3=123321, # 4 = 1 234432 1, ... 

By adding to the all elements of the analytical formula 1(1) from the right 
the operator A ~ k , having the same form as one located from the left, we obtain 
the new formula of the Magic square 3x3. This formula allows easy to construct 
examples of Magic squares 3x3 both from numbers of the first subsequence {see 
Fig. 2(1)} and from numbers of the second subsequence {see Fig. 2(2)}. 



171819191817 


101 1 12121 1 10 


151617171615 


121314141312 


141516161514 


161718181716 


131415151413 


181920201918 


111213131211 



( 1 ) 



17181920191817 


10111213121110 


15161718171615 


12131415141312 


14151617161514 


16171819181716 


13141516151413 


18192021201918 


11121314131211 



( 2 ) 



Fig. 2. Constructing Magic squares 3x3 from ^-truncated numbers of S,-series. 
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3. Magic squares 3x3 in size from k-truncated numbers of Sy-series. By 
comparing numbers of S 3 -series {see (11)} and ^-series {see point 2} with each 
to other one can conclude that numbers of ^-series resemble numbers of the 
first subsequence of ^-series and distinguish from them on the order of the 
natural numbers movement. The example of the Magic square 3x3 from 
numbers of Sj-series is presented in Fig. 3. This square is constructed by means 
of methods described in point 1 and 2. Thus, in spite of the mentioned 
difference between numbers of Sj-series and S 2 -series, the methods, discussed 
above, can be applied for solving problems on constructing Magic square 3x3 
from numbers of iSj-series. 



201918181920 


131211111213 


181716161718 


151413131415 


171615151617 


191817171819 


161514141516 


212019192021 


141312121314 



Fig. 3. Constructing Magic squares 3x3 from k-truncated numbers of S 3 -series. 



4. Magic squares 3x3 in size from k-truncated numbers of Srseries. In 
contrast to considered Smarandache sequences the digit 0 is absent in numbers 
of S 4 - series. Besides, the order of the movement for digits 1, 2, .... 9 can not be 
changed and after digit 9 can be the only digit 1. These peculiarities of numbers 
of S 4 - series make too difficult the solving problems on constructing Magic 

square 3x3. It is evident that by using A~ k -operator one can easy construct 
classical square 4(1) {the Magic square of natural numbers from 1 to 9}. Since 
by means of A ~ k -operator such square can be constructed from numbers of 
any Smarandache sequence {for instance, see (6)} , the example of the square 
4(1) is banal. The example of the Magic square 3x3, presented in Fig. 4(2, 3), is 

less trivial. 



78 


1 


56 


23 


45 


67 


34 


89 


12 



a A' 2 

4 




A" 1 *, 


*2 


a j A” 1 


A -2 a 7 A -3 


A' 1 a 5 A -2 


A -1 a A -1 

4 


A -3 a 6 A -1 



8 


1 


6 


3 


5 


7 


4 


9 


2 



(1) (2) (3) 

Fig. 4. Constructing Magic squares 3x3 from k-truncated numbers of S 4 -series. 
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5. Magic squares 3x3 in size from k-truncated numbers of S 5 -series. As 
compared with another Smarandache sequences of the 1st kind the numbers of 
*S 5 -series {see (16)} have the following peculiarity: the circular permutation of 

natural numbers is allowed in them. The analytical formula of Magic square 
3x3, presented in Fig. 5(1), is just constructed with taking into account the 
pointed peculiarity of discussed numbers. Examples of the Magic square 3x3, 
obtained from formula 5(1) at n = 2, 3 and 4, are depicted in Fig. 5(2, 3, 4) 
correspondingly. By analysing these squares it is easy to find more simple form 
of the analytical formula 5(1) {see Fig. 5(5), where a r-1 is the ( n -l)th number 

of <S r series, M is general amount of digits in the number a }. 



A 1<2 n(n+ 15)/2 + 21 


a n(n+\)/2 


A a /»(n+ll)/2 + 10 


_2 

A a n(n+5)/2 + 1 


A n 

a n(n+9)/2 + 6 


A a n(n+l3)/2 + 15 


A 3fl n(n+ 7)/2 + 3 


A n 

°/r(/z+17)/2 + 28 


A la n(n+ 3)12 



( 1 ) 



11123 


4123 


9123 


6123 


8123 


10123 


7123 


msm 


5123 



1012 


312 


812 


512 


712 


912 


612 


1112 


412 



91 


21 


71 


41 


61 


81 


51 


101 


□ 



( 2 ) 



( 3 ) 



(4) 



(*+ 7) 10 M+ a n _y 


n\b M +a , 

ft— 1 


(n+5)l0 M +a n _ l 


(n+2)10 A/ +a n _, 


(n+4)lQ M + a n _ { 


(n+6)I0 M +a n _, 


(*+3)10 M +a n _ x 


(/i+8)10 M +a n _i 


{n+\)\<3 M + a n _ y 



( 5 ) 



Fig. 5. Constructing Magic squares 3x3 from ^-truncated numbers of S 5 -series. 
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6. Magic squares 5x5 in size from k- truncated numbers of S ^series. Numbers 
of 5^-series {see (18)} resemble both numbers of the first subsequence of S 2 - 
series and numbers of S 3 -series {see points 2 and 3} .The example of the Magic 
square 3x3 from numbers of S 6 -series is presented in Fig. 6. This square is 
constructed by means of methods described in points l - 3. Thus, in spite of 
the mentioned difference between numbers of «S 6 -series and S 2 - y ^-series for 
solving problems on constructing Magic square 3x3 from numbers of S 6 -series 
the methods, discussed above, can be applied. 



2527293132302826 


1 1 13151718161412 


2123252728262422 


1517192122201816 


1921232526242220 


2325272930282624 


1719212324222018 


2729313334323028 


1315171920181614 



Fig. 6. Constructing Magic squares 3x3 from k - truncated numbers of S 6 -series. 
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Perfect Powers in Smarandache Type Expressions 



Florian Luca 

In [2] and [3] the authors ask how many primes are of the Smarandache 
form (see (10]) x y + y x , where gcd ( x , y) = 1 and x, y > 2. In [6] the author 
showed that there are only finitely many numbers of the above form which 
are products of factorials. 

In this article we propose the following 

Conjecture 1. Let a , b, and c be three integers with ab £ 0. Then the 
equation 

ax y + by x = cz 11 with x, y, n> 2, and gcd (x, y) = 1, (1) 

has finitely many solutions (x, y, z, n). 

We announce the following result: 

Theorem 1. The ”abc Conjecture” implies Conjecture 1. 

The proof of Theorem 1 is based on an idea of Lang (see [5]). 

For any integer k let P(k) be the largest prime number dividing k with 
the convention that P(0) = P(±l) = 1. We have the following result. 

Theorem 2. Let a, b, and c be three integers with ab ^ 0. Let P > 0 
be a fixed positive integer. Then the equation 

ax y +by x = cz n with x, y, n> 2, gcd (x, y) = 1, and P(y) < P, (2) 

has finitely many solutions (x, y, z, n). Moreover, there exists a com- 
putable positive number C depending only on a, b, c, and P such that all 
the solutions of equation (2) satisfy max (x, y) <C. 

The proof of theorem 2 uses lower bounds for linear forms in logarithms 
of algebraic numbers. 

Conjecture 2. The only solutions of the equation 

x y ±y x = z n with x, y, n > 2, z > 0, gcd (x, y) = 1, (3) 

are (x, y, z , n) = (3, 2, 1, n). 

We have the following results: 

Theorem 3. The equation 

x y ±y x = z 2 with x, y > 2, and gcd (x, y) = l, (4) 

has finitely many solutions (x, y, z) with2 \ xy. Moreover, all such solutions 
satisfy max (x, y) < 3 • 10 143 . 




The proof of Theorem 3 uses lower bounds for linear forms in logarithms 
of algebraic numbers. 

Theorem 4. The equation 

2 y + y *=z n (5) 

has no solutions ( y , z, n } such that y is odd and n > 1. 

The proof of theorem 4 is elementary and uses the fact that Z[i>/2] is 
an UFD. 

2. Preliminary Results 



We begin by stating the abc Conjecture as it appears in [5]. Let A; be a 
nonzero integer. Define the radical of k to be 

No(k) = n P (6) 

p 1 fc 

i.e. the product of the distinct primes dividing k. Notice that if x and y are 
integers, then 

No(xy) < No(x)No(y), 
and if gcd ( x , y) = 1, then 

No{xy) = No(x)N 0 (y). 

The abc Conjecture ([5]). Given e > 0 there exists a number C(e) 
having the following property. For any nonzero relatively prime integers 
a , b, c such that a + b = c we have 

max(|a|, \b\, |c|) < C(e)N 0 (abc) l+€ . 

The proofs of theorems 2 and 3 use estimations of linear forms in logarithms 
of algebraic numbers. 

Suppose that Ci, — , 0 are algebraic numbers, not 0 or 1, of heights not 
exceeding Ai, ..., Ai, respectively. We assume A m > e e for m = 1, ..., 1. 
Put Q = log Ai~.log .4*. Let F = Q(Ci, ..., 0]- Let •••> be integers, 
not all 0, and let B > max In^. We assume B > e 2 . The following result 
is due to Baker and Wiistholz. 

Theorem BW ([!]). Iftf'-C? 1 £ 1, then 

1C 1 -C ~ 1| > \ exp(— (16(Z + l)d F ) 2 ( ,+3 )D log B). (7) 

In fact, Baker and Wurtholz showed that if log Ci, —, logCz are any 
fixed values of the logarithms, and A = n\ log £i + ... + n/ log Q -fc 0, then 

log |A| > -(16 Zd F ) 2(f+2) D log B. (8) 
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Now (7) follows easily from (8) via an argument similar to the one used by 
Shorey et al. in their paper [9]. 

We also need the following p-adic analogue of theorem BW which is due 
to Alf van der Poorten. 

Theorem vdP ([7]). Let n be a prime ideal of F lying above a prime 
integer p. Then, 



ord.(C-Cr -1) < (l6(l + l)<iF) 12<1+1) i^n(IogB) 2 . (9) 

lu O r 



We also need the following two results. 

Theorem K ([4]). Let A and B be nonzero rational integers. Let 
m >2 and n > 2 with mn >6 be rational integers. For any two integers x 
and y let X = max (|x|, |y|). Then 

P(Ax m + By n ) > C (log 2 X Iog 3 X) 1/2 (10) 

where C > 0 is a computable constant depending only on A, B, m and n. 

Theorem S ([8]). Let n > 1 and A, B be nonzero integers. For integers 
m > 3, x and y with |x| > 1, gcd (x, y) = 1, and Ax m 4- By n ^ 0, we have 

P(Ax m + By n ) > C((log m)(log log m)) 1/2 (11) 

and 

| Ax m + By n | > exp^C ((log m) (log log m)) 1/2 j (12) 

where C > 0 is a computable number depending only on A, B and n. 

Let K be a finite extension of Q of degree d, and let O k be the ring 
of algebraic integers inside K. For any element 7 € Ok, let [7I be the ideal 
generated by 7 in Ok- For any ideal / in Ok, let N(I) be the norm of I. 
Let 7Ti, 7 t 2> be a set of prime ideals in Ok- Put 

p = max P(lV(ir i )). 

Write 

= N for i •••> 1 

where pi, P2, .», pi € Ok and h is the class number of K. Denote by <S the 
set of all elements a of O k such that [a] is exclusively composed of prime 
ideals 7 Ti, ^2, ttj. Then we have 
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Lemma T. ([9]). LetaeS. Assume that 

[a] = Tr b 1 1 ir b 2 \..ir b l ‘. 

There exist a 0 € Ok with \N(0)\ < , p dhl and a unit e e Ok such that 

a = epvTvV -VT ■ 



Moreover, 



bi = aih + Ci for some 0 < d < h. 



3. The Proofs 

The Proof of Theorem 1. We may assume that gcd (a, b, c) = 1. 

By Ci, C 2 , we shall denote computable positive numbers depending 

only on a, b, c. Let (x, y, z, n) be a solution of (1). Assume that x > y, 

and that x > 3. Let d = gcd (ax y , by x ). Notice that d \ ab. Equation (1) 

becomes „ , _ „ 

ax y by _ cz 

~T + ~d ~ ~d' 

By the abc Conjecture for e = 2/3 it follows that 

max (|ox I '|, \by x \ \cz n \ ) < C ( 2 / 3 )^( a&c ) . N ^ xyz )5/3 

Let 

Ci = C(2/3)N Q {abc) 5/3 

Since d > 1 , and |6[ > 1 , from inequality (14) it follows that 
y x < \by x \ < Ci(xy\z\f/ 3 < Cix i0 l 3 \z\V z . 

Since x > min (y, 3), it follows easily that y x ■> x y . Hence, 



W n = 



-x y + -y x < Ciy* 
c c 



where C 2 = ^ . We conclude that 

\c\ 

\z\ < C\ ,n y x/n < C\ ,2 y x/n . 
Combining inequalities (15) and (16) it follows that 

y* < C,Cf / V 0/ ¥ 5 * /3 " ) , 

or 

y x(l—5/Zn) < C 3 x 10/ 3 , 



(13) 



(14) 



(15) 



(16) 



(17) 
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( 18 ) 



where C 3 = CiCf /6 . Since 2 < y and 2 < n, it follows that 

2 i /6 < 2 x(l-5/3n) < q^IO/Z 



Inequality (18) clearly shows that x < C 4 . 

The Proof of Theorem 2. We may assume that 



P > max (P(a), P(6), P(c)). 



By Ci, Ci, we shall denote computable positive numbers depending 
only on a, b, c, P. We begin by showing that n is bounded. Fix d € 
{2, 3, P - 1}. Suppose that x, y, z, n is a solution of (2) with n > 3 
and d | y. Since 

bif = cz n - a(x y / d ) d (19) 

it follows, by Theorem S, that 

P = P(by x ) = P(cz n - a(x y/d ) d ^j > C x ((logn)(log logn)) 1/2 (20) 

where C x is a computable number depending only on a, c, d. Inequality 
(20) shows that n < Ci- 

Suppose now that ny > 6. Let X = max (x, |z|). From equation (19) 
and theorem K, it follows that 

P = P(by *) = P(cz n — ax y ) > C 3 (log 2 X log 3 X) 1/2 , (21) 

where C 3 > 0 is a computable constant depending only on a, c, and Ci- 
From inequality (21) it follows that X < C 3 . Let C\ = max (C2, C 3 )- It 
follows that, if ny > 6, then max (x, \z\, n) < C*. We now show that y is 
bounded as well. Suppose that y > max (C4, e 2 ). Rewrite equation (2) as 



M n _ 

lalx 5 




( 22 ) 



Let A > e e be an upper bound for the height of -b/a and C 4 . Let Q = 
(log A) 3 . From theorem BW we conclude that 



log |c| + n log \z\ - log |a| - y log x > - log 2 - 64 12 ft log y. (23) 



Since x > 2, and max (x, \z\, n) < C 4 , it follows, by inequality (23), that 



y log2-64 12 Qlogy < y log x-64 12 Q logy < C 4 logC 4 -log|a|+log|c|+log2. 

(24) 

From equation (24) it follows that y < C 5 . 
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Suppose now that n = y = 2. We first bound z in terms of x. Rewrite 
equation (2) as 

l&l 



z 2 = 2 X 



1 + 



(iXIOI- 



(25) 



Let Cg > 0 be a computable positive number depending only on a and b 
such that 



a /x \ i 
6 v2*/ < 2 



for x > Cg. 



(26) 



From equation (25) and inequality (26), it follows that 

r*2 ^ ^ _ |^l/ l ft l / ft»2 



^l<^l( 1 H»<^< 2 m i+ lfl(F))< 2 11l <-> 

for x>Cg. Taking logarithms in inequality (27) we obtain 

xCr + Cg < log z < xCr + Cg for x > Cg (28) 

where C 7 = !?§2, C a = lo 8 i&! ~ lo S 2M , ^ c , = W ~ '°S N . We 
now rewrite equation (2) as 



(cz) 2 — acx 2 = ab2 x . 



(29) 



Let a = y/ac. Then 



{cz + ax){cz — ax) = cb2 x . (30) 



We distinguish 2 cases. 

CASE 1. ac < 0. Let K = Q[a]. Since ac < 0, it follows that all the 
units of Ok are roots of unity. Since K is a quadratic field, it follows that 
the ideal [2] has at most two prime divisors. Since 

gcd {[cz + ax], [cz — ax| j | 2[abc] 

it follows, by lemma T, that 

cz + ax = e/?p“ (31) 

where ^ - 1 < u < and e,0, p € O k are such that |e| = 1, |p| = 2 h/ 2 , 

where h is the class number of K, and |/3| < Cio where Cio is a computable 
number depending only on a, b, and c. Conjugating equation (31) we get 

cz — ax = epp*- (32) 

From equations (31) and (32) it follows that 

2qx = ePp u { 1 - (-e- 2 )(/?r 1 £(p)-“(p) u ). 
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Hence, 



(33) 



2M* = i/3ii P r|i - {-e- 2 m- x m~ u {vY | 

Taking logarithms in equation (33) we obtain 

log(2|a|)+logx = log|/?| + ulogp + log|l-(-e -2 )(/?) _1 /?(p)"“(p)“|. (34) 

Let A, and P be upper bounds for the heights of —e~ 2 (P)~ l P and p, respec- 
tively. Assume that min (A, P) > e e . Let Q = log A(logP) 2 . Assume also 
that — > 1 + e 2 . From equation (34), theorem BW, the fact that |p| = 2 h / 2 , 

h x x 

and the fact that --l<u<-, we obtain that 
h h 

log(2|a|) + log x > log \P\ + u log |p| - log 2 - 64 12 Q log u > 

log \p\ + - l) • (j) log 2 - log2 - 64 12 0 log (x/h). (35) 

Inequality (35) clearly shows that x < Cn- 

CASE 2. ac > 0. We may assume that both a and c are positive. If 
6 < 0, equation (2) can be rewritten as 

|a|x 2 - \b\2 x = \c\z 2 > 0 (36) 

Equation (36) clearly shows that x < Cn- Hence, we assume that 6 > 0. 
We distinguish two subcases. 

CASE 2.1. y/ac € Z. In this case, from equation 

(c\z\ + ox){c\z\ — ox) — bd2 x 



and from the fact that 
gcd 

it follows easily that 




+ ax, c\z\ — axj \ 2 acb 


(37) 


c\z\ + ax = p2 u 
c\z\ — ax = 7 


(38) 



where /?, 7, u are positive integers with 0 < P < be, 7 < (be) ■ (2acb) and 
u > x — ord2(2ac6). From equation (38) it follows that 



2o:x = p2 x - 7. (39) 

From equation (39), and from the fact that 0 < p < be, 7 < (be) ■ (2 cceb), 
and u > x — ord2(2ac6), it follows that x < C13. 



69 




CASE 2.2. \fac £ Z. Let K = Q[a]. Let e be a generator of the torsion 
free subgroup of the units group of Ok- From equation (37) and lemma T, 
it follows that 

c|z| + ocx — e m 0iPi (40) 

X X 

where - — 1 < u < — , and 0, pi € Ok are such that 1 < 0i < Cu for some 
h h 

computable constant Cu, and 1 < p\ < 2 -e. From equation (40), it follows 
that 

c|z| -ax = e~ m 02P2 (41) 

where 02 = \0i\ 2 /0i, and p 2 = 2 A /Pi- Suppose now that x > Ce- Since 

e m = 1 (c|z| + ax) 

X X 

it follows, from inequality (28), and from the fact that — — 1 < u < — and 

h h, 

1 < pi < 2 h ■ e, that 

\m\ < Cisx + Cie for x > C$, (42) 



for some computable constants C is and C\e depending only on a, b, and c. 
From equations (40) and (41), it follows that 

2 ax = e m 0iPi ■ (l - e -2m (/?i ) “ 1 /?2 (p i ) ~“P 2 ) 

or / \ 

2ax = (c|z| + ax) • (l - e- 2m {0{)- 1 0 2 {pi . (43) 

Let Ai, A 2 , A3, A 4 be upper bounds for the heights of e, (0i)~ 1 02, Pi, P2 
respectively. Assume that min (Ai, A2, A3, A4) > e®. Denote Cl = 
ntilogA,. Denote Cu = max (2 Ci 5 , l/h). From inequality (42), it 
follows that 

max (2|m|, u) < Cux + Cie- (44) 

Let B = Cux + Ci6- Taking logarithms in equation (43), and applying 
theorem BW, we obtain 



log(2a) + logx = log(c|z| 



1 j 






log(c|z| + ax) - log 2 - 80 14 fi log(C 17 x + C 16 ). (45) 

Combining inequalities (28) and (45) we obtain 

log(4a) +logx + 80 14 fHog(Ci7X + C 16 ) > log(c|z| +ax) > logz > C 7 x + Cs 
This last inequality clearly shows that x < Cis- 
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The Proof of Theorem 3. We treat only the equation 



x y + y x = z 2 . 



We may assume that x is even. First notice that, since gcd ( x , y) = 1, it 
follows that gcd (x, z) = gcd ( y , z) = 1. Rewrite equation (4) as 

X" = ( 2 + S' /J )(^-S I/2 ). 



Since gcd (z, y x / 2 ) = 1 and both z and y are odd, it follows that 

gcd (z + y l/2 , z — y x/2 ) =2. 

Write x = 2d x di such that either one of the following holds 

f z + y x/2 = 2 y_1 d^ Qr f z + y x/2 = 2djf 

\ z _ y */2 = 2 d 2 ° r \ z - y I / 2 = 2“- 1 4 



Hence, either 

y 1 / 2 = 2^- 2 4 - 



(46) 

(47) 



y 1 / 2 = d\- 2 y ~ 2 <% (48) 

We proceed in several steps. 

Step 1. (1) If x > y then either y < 9 and x < 27, or y > 9 and 
x < 3 y. 

(2) //x < y and y > 2.6 • 10 21 , then y < 4x. 

(1) Assume first that x > y. Since 

y x/2 = 2 y-2 d y _ ^ or y x/2 = rf y _ 2 y-2^ 



it follows that 

y s /2 < 2 y ~ x d\ < (2di) y < x y or y x ^ 2 < <£{ <x y . 

Hence, 

-logy < ylogx. 

Inequality (50) is equivalent to 

x y 

<2 

log x logy 



(49) 

(50) 



(51) 



If y < 9, then one can check easily that (51) implies x < 27. Suppose now 
that y > 9. We show that inequality (51) implies x < 3y. Indeed, assume 
that x > 3y. Then 



3y _ 3y < _x_ < _2y_ 
log 3 + logy log(3y) — logx logy 



(52) 
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Inequality (52) is equivalent to 

3 log y < log 9 + 2 log y 

or y < 9. This contradiction shows that x < 3y for y > 9. 
(2) Assume now that x < y. Suppose first that 

y x/2 = 2 y-2 d y _ ^ 



In this case 



(2d{) y >2y- 2 dX=<% + y x ' 2 ><% 
Sin< 

Suppose now that 



therefore 2di > d 2 . Since x = 2did 2 , it follows that 2di > y/x, or di > 






In this case, 

d\ > 2 y - 2 <% > 4 

or d\ > d 2 . We obtain that di > y/did 2 = 

If equality (47) holds, it follows that 

y x/2 = 2 y- 2 dX\i -2 _(y_2) (^-) y | > d?|i - 2_(y_2) (^) y |- ( 53 ) 

On the other hand, if equality (48) holds, then 

S*/ 2 =<|l-2“- 2 (^) , '|. (54) 



From inequality (53) and equation (54), we conclude that, in. either case, 



v zn > <9 






(55) 



for some e € {±1}- Suppose now that x > e e . By theorem BW, and 
inequality (55), it follows that 

^ logy > ylogdt - log 2 - 48 10 e log x log y > 



or 



yfx 



y log log 2 - 48 lu e log x log y 



in X , , y/x 

48 10 e log x log y + log 2 + - logy > ylog — . 



(56) 

(57) 
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CASE 1. Assume that x < 2 6 . From inequality (57), it follows that 

f~e. 

48 10 e • 61og2 • logy 4- log 2 4- 2 s logy > ylog | 



or 



or 



(48 10 e • 6 log 2 4- 2 s ) log y + log 2 > | 

2(48 10 e • 6 log 2 + 2 5 + 1) > T — . 
v logy 



(58) 



Let Ci = 2(48 10 e • 6 log 2 + 2 s 4- 1). From inequality (58) and lemma 2 in 
[6], it follows that 



y < Ci log 2 Ci < 2(48 10 e - 6 Iog2 + 2 5 + 1) • 42 2 < 2.6 • 10 21 . 
CASE 2. Assume that x > 2 6 . Then, 

*>^2 Vi. 

Inequality (56) becomes 
>xo 



(59) 



or 



or 



or 



x 1 

48 10 e log x log y + log2 + - logy > - ylogx 



3e48 10 log x log y + log 8 + - xlogy > ylogx 



(3e48 10 4- 1) log x logy + - x logy > y logx 
3e48 10 + 1 + 1 t ~ — > V 



2 logx logy' 



Assume first that 



In this case, 



3 x 



< 3e48 10 + 1. 



2 logx 



i^<f ( 3e48l0 + 1 )' 



(60) 

(61) 

(62) 



Let C-i — - (3e48 10 + l). From inequality (62) and lemma 2 in [6], it follows 



that 



x < C 2 log 2 C 2 < ^(3e48 10 + 1) - 41 2 < 6 • 10 20 . 

O 

In this case, from inequalities (60) and (61), it follows that 



< 2(3e48 10 + 1). 



logy 



(63) 



(64) 
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Let Cz — 2(3e48 10 + 1). It follows, by inequality (64) and lemma 2 in [6], 
that 

y < C 3 log 2 Cz < 2(3e48 10 + 1) • 42 2 < 1.8 • 10 21 . (65) 

Assume now that y > 2.6 -10 21 . From inequality (59), it follows that x > 2 6 . 
Moreover, since inequality (65) is a consequence of inequality (61), it follows 
that 

- — > 3e48 10 + 1. (66) 

2 log x ~ 

From inequalitites (60) and (66) it follows that 



y 

log x log y' 



(67) 



We now show that inequality (67) implies y < Ax. Indeed, assume that 
y > 4a;. Then inequality (67) implies 

3a: y Ax _ Ax 

logx logy “ log(4x) log x + log 4 



or 

3 log x + 3 log 4 > 4 log x 

or 3 log 4 > logx which contradicts the fact that x > 2 6 . 
Step 2. Ify>3 ■ 10 143 , then y is prime. 

Let 



yX/2 = 2 y-2 d y _ jy or y-/2 = d y- 2 v- 2 4_ (68) 



Notice that if y 1 / 2 = 2 y ~ 2 d\ - d%, then gcd (2d!, d 2 ) = 1. Let p | y be a 
prime number. Since p / 2did 2 = x, it follows, by theorem vdP, that 

| < max (ord p (2«- 2 d?-d|), ord p (d?-2^" 2 4)) < 48 36 e^ log 2 ylogx. 
2 ' ' ° g (69) 

By step 1, it follows that 

-y < x < 2 • 48 36 e— log 2 y log(4y) < 4 • 48 36 e— log 3 y. (70) 
4 ~ logp logp 



Hence, 

-4- < 16 • 48 36 er^- < 16 • 48 36 ep. (71) 

log 3 y logp 

Suppose that y is not prime. Let p | y be a prime such that p < yfy. From 
inequality (71) it follows that 



< 16 • 48 36 e 

log y 
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or 



Vv 



log (y/y) 



< 128 • 48 36 e. 



(72) 



Let k = y/y and C 4 = 128 • 48 36 e. By inequality (72) and lemma 2 in [6], it 
follows that 



y/y = k < C 4 log 4 Ci = 128 • 48 36 e • 146 4 < 5.3 • 10 



^1 



or 



y < (5.3 • 10 71 ) 2 < 3 • 10 1 
This last inequality contradicts the assumption that y > 3- 10 143 . 

Step 3. If y > 3 • 10 143 , then x>y. 

Let y = p be a prime. If y x/2 = 2 y ~ 2 d\ — d%, it follows, by Fermat’s 
little theorem that 

2~ 1 di -d 2 = 2 y ~ 2 d\ -d%= y x/2 = 0 (mod p), 

therefore 



■>143 



(73) 

(74) 



di = 2d 2 (mod p). 

On the other hand, if y 1 / 2 = d\ — 2 y ~ 2 d%, then 

di - 2 -1 d 2 = d\ — 2 y ~ 2 d\ = y x/2 = 0 mod p), 

therefore 

d 2 = 2d\ (mod p). 



(75) 



(76) 



Suppose that x < y. From congruences (75) and (76), we conclude that, in 
both cases, x is a perfect square. Hence, 



»’ = z 2 - (v'J) 2 '' = (* + (v®’) ■ {z - (,/i)')- 
From equation (77) it follows that 

(z-(VZ) y = 1 

\ z + (>/x) 9 = y x 

Hence, 

2{y/x) y = y x - 1. 

It follows, by equation (79) and theorem BW, that 



(77) 



(78) 



(79) 



0 = log 



y x - 2{y/x) v = log(y x ) + log 1 — 2y x (/x) 1 
xlogy- log 2 - 64 12 e log 2 y log x. 



> 



( 80 ) 
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From inequality (80) and Step 1 it follows that 

log 2 + 64 12 e log 3 y > x log y > - - 

or 

4 log 2 + 4 • 64 12 e log 3 y > ylogy 
or 

(4 • 64 12 e + l) log 2 y > y 
or 

4-64 12 e + l > — (81) 
log y 

Let Cs =4- 64 12 e + 1. By inequality (81) and lemma 2 in [6j it follows that 
y < C 5 log 3 C 5 < (4 • 64 12 e + 1) • 53 3 < 8 • 10 27 . (82) 

The last inequality contradicts the fact that y > 3 • 10 143 . 

Step 4. Suppose that y > 3 • 10 143 . Let y = p be a prime. Then, 
with the notations of Step 1, every solution of equation (4) is of one of the 
following forms: 

(1) y x/2 = 2 v ~ 2 d\ — d .2 with y = p, d x = 2+p, d 2 = 1, x = 4 + 2p 

(2) y x/2 — d\ — 2 y ~ 2 d% with y — p, d\ = d 2 = 1, x = 3p - 1 

(3) y x/2 =d\- 2 y ~ 2 d% with y = p, d x = - - 1 , d 2 = 3, x = 3p-9 

We assume that y > 3 • 10 143 . In this case, y = p is prime, and x > y. 
From Step 1 we conclude that x < 3 y. Moreover, from the arguments used 

Jx 

at Step 1 it follows that d x > Since x = 2did 2 , it follows that 

d 2 < y/x < y/Z y = y/Zp. 

By the arguments used at Step 3 we may assume that x is not a perfect 
square. We distinguish the following cases. 

CASE 1. d 2 — 1. By congruences (75) and (76) it follows that d\ = 
2 (mod p), or 2<2i = 1 (mod p). 

Assume that d x = 2 (mod p). Since x = 2d x , and p = y < x < 3y = 3p, 
it follows that d\ = 2+ p and x --- 2d x =4 + 2 p. 

Assume that 2d x = 1 (mod p). Again, since x = 2di, and p = y <x < 

3 y = 3 p, it follows that d\ = — , and i = 3p - 1. 

CASE 2. d 2 = 2. By congruences (75) and (76) it follows that d\ = 
4 (mod p), or di = 1 (mod p). One can easily check that there is no solution 
in this case. Indeed, if d x = 4 (mod p), it follows that d x > p + 4. Hence, 
x = 2d\d 2 > 4(p + 4) > 3p = 3y which contradicts the fact that x < 3y. 
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Similar arguments can be used to show that there is no solution for which 
d 2 = 2 and d x = 1 (mod p). 

CASE 3. d 2 = 3. By congruences (75) and (76) it follows that di = 
6 (mod p), or 2 d x = 3 (mod p). One can easily check that there is no 

solution for which di = 6 (mod p). Suppose that 2d\ = 3 (mod p). Since 

p — 3 

p = y<x<3y = 3p and x = 2did 2 = 6di, it follows easily that di = , 

and x = 3p — 9. 

CASE 4. d 2 = fc > 4. 

If it is even, then, by congruences (75) and (76), it follows that d x = 
2k (mod p), or d x = k/2 (mod p). Since x is not a perfect square it follows 
that d x > p+k/2, therefore x > 2pk+k 2 > pk > 4p > 3p = 3y contradicting 
the fact that x < 3 y. 

If k is odd, then, by congruences (75) and (76), it follows that di = 

rQ jjj* 

2k (mod p), or 2d x = k (mod p). We conclude that d x > — - — , therefore 

x = 2did 2 > k(p—k). Since k{p-k ) > 3 p for 5 < k < ,/3p and p > 3- 10 143 , 
we conclude that x> 3p = 3y contradicting again the fact that x < 3 y. 

Step 5. There are no solutions of equation (2) with y > 3 - 10 142 and 
x even. 

According to Step 4 we need to treat the following cases. 

CASE 1. 

yi/2 _ 2 y ~ 2 d\ — d| with y = p, di = 2 + p, d 2 = 1, x = 4 + 2 p. (83) 

Hence, 

p2+ P = 2 p- 2 (2 + p) p - 1 > 2 p-3 (2 + p) p . (84) 

Taking logarithms in inequality (84) we obtain 

(2 + p) log p > (p - 3) log 2 + p log(p + 2) 



or 

2 log p + p(logp - log(p + 2)) > (p - 3) log 2. 
It follows, by inequality (85), that 



2 log p > (p — 3) log 2 



(85) 



or 

plog2 < 21ogp + 31og2 < 5 logp. (86) 

Inequality (86) is certainly false for p = y > 3 • 10 143 . 

CASE 2. 

3p — 1 

y*/ 2 = d y - 2 y_2 df with y = p, di = t> ~ , d 2 = 1, x = 3p - 1. 
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Hence, 



0r O 

P (p_1)/2 < (f ) P - (87) 

Taking logarithms in inequality (87) it follows that 

p — 1 

— — logp<plogl.5 

or 9 

logp < 7? - log 1.5 < 3 log 1.5 < log 1.5 3 . 

p- 1 

It follows that p < 1.5 3 < 4 which contradicts the fact that p > 3 • 10 143 . 
CASE 3. 

y l/2 = d\ - 2 y ~ 2 d.2 with y = p, di = - - \ d 2 = 3, x = 3p - 9. 
Hence, 

p(3 P -9)/2 = - 2 P-2 3 P < ? < pP- (88) 

3p — 9 

From inequality (88) it follows that — - — < p or p < 9 which contradicts 
the fact that p = y > 3 • 10 143 . 

The Proof of Theorem 4. The given equation has no solution 
(y, z, n) with n > 1 and y odd, y < 5. Assume now that y > 5. We 
may assume that n is prime. We first show that n is odd. Indeed, assume 
that (y, z) is a positive solution of y 2 + 2 y = z 2 with both y and z odd. 
Then (z + y)(z - y) = 2 y . Since gcd (z + y, z - y) = 2 it follows that 
z —y = 2 and z + y = 2 y_1 . Hence, y = 2 y-2 — 1. However, one can easily 
check that 2 y-2 — 1 > y for y > 5. 

Assume now that n = p > 3 is an odd prime. Write 

(y + 2 (y-1)/2 • i'/2) • (y - 2 (y_1)/2 • iy/Tj = z n 

Since Z[tV2] is euclidian and 

gcd (y + 2 (y ~ 1)/2 • i'J 2, y - 2 (y - 15 / 2 • 1 V 2 ) = 1 
it follows that there exists a, b € Z such that 

f y + 2 (y-1) / 2 - i\/2 = (a + bi-v/if) 

; L (89) 

1 y — 2( y-1 ^ 2 • i^/2 = — biV2j 
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